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1 Introduction 

1.1 Josephson junction arrays 

Two-dimensional (2D) Josephson junction arrays {J J A) offer a unique opportunity for study- 
ing a variety of topics in 2D physics, such as phase transitions, non-linear dynamics, perco- 
lation, frustration and disorder, in relatively <Cclean^ experimental realisation. Fabrication 
of arrays and their basic physical properties have been described in various articles [1], [2]. 
A JJA consists of islands, positioned in some periodic or irregular arrangement, that be- 
come superconducting below a given transition temperature T°. Below this temperature 
each island I is characterized by its superconducting wave function 

i^ = me''' (1.1) 

with its amplitude and its phase. For all matters and purposes one can assume that has 
the same value in each island, such that the phase is the only relevant variable. The islands 
are linked to each other by the Josephson coupling. The potential energy of the array is then 
given by 

H = Y,Ml-cos{ei-0i,)] (1.2) 

(lU) 



The sum can usually be restricted to nearest neighbors in the array and the Josephson 
coupling J is related to the critical current Ic by 

./ = |:4 (1.3) 

If the array is placed in a perpendicular magnetic field the phase differences in (1.2) have to 
be replaced by the <C gauge invariant ^ combination 

4>w =0i- Ov - Alv (1.4) 

where Aui is the line integral 

^11' = ^ I dr-A(r) (1.5) 

of the applied vector potential. Most of our subsequent calculations will be based on the 
assumption of a large magnetic penetration depth. Then any vector potential appearing in 
our equations describes an applied external electromagnetic field, whereas internal field fiuc- 
tuations are neglected. 



When the phases of neighboring islands differ from each other a supercurrent 

Kiif = Jill sin0/// (1.6) 
is flowing through the junction (/, I') and the two islands get charged. The electrostatic energy 

^ w 

has to be added to (1.2). Here Qi is the charge of island related to the excess number 5Ni 
of Cooper pairs by Qi = 2e6Ni , and Cui is the capacitance matrix of the array. Its diagonal 
elements determine the charging energy with respect to ground and the nearest neighbor 
elements represent the junction capacitance. The excess number 6Ni can be identified with 
a canonical momentum Pi conjugate to the phase [1,3,4]: 

[Pi,ei,] = -6ii, (1.8) 

The electrostatic energy (1.7) thus plays the role of a (non-local) kinetic energy that has to 
be added to the potential energy (1.2) yielding thus the complete Hamiltonian 

H = ^J2('^efPiiC-%>Pi, + Y.Jii,{l-cos(l)ii,) (1.9) 

W {W) 

of the array on which our equations of motion will be based in section 2. 



Much work has been done in order to elucidate the thermodynamic properties of (1.9) 
which is indeed a challenging problem of statistical mechanics. The simplest case taking into 
account only the potential energy (1.2) in zero applied field is the classical XY-model in 2 
dimensions. In its ground state all the <Cspins^> are parallel, respectively all the phases Oi 
are equal. This system allows to study the critical behaviour of a 2D system with a com- 
plex order parameter. It gives rise to the well-known Berezinski-Kosterlitz-Thouless phase 



transition [5,6] bringing in the basic concept of vortex excitations. The basic features are 
summarized in the following subsection. 



The full form (1.9) adds two complications: 

a) A J J A in an applied magnetic field corresponds to the <Cfrustrated^ XY-model. 
Owing to the additional term Aui in the coupling energy the phase configuration in the 
ground state is more complicated - corresponding to magnetic field induced vortices - and 
the energy per bond is less negative than —2 J. The effect of Aui is measured in terms of the 
flux threading through an elementary cell of the array, the <Cfrustration parameter^ / 



the sum running over the bonds going around the cell. For simple - rational - values of / the 
ground state phase pattern is known (see [1] and the references given there). For certain /- 
values interesting infinite degeneracies have been found in triangular arrays [7] The behavior 
at finite temperatures, in particular the various kinds of phase transitions that a frustrated 
array is supposed to undergo is less well understood [1]. 

b) Owing to the commutation relations (1.8) the full Hamiltonian (1.9) describes a quan- 
tum system, sometimes called <C quantum XY-model^ although its operator algebra is, of 
course, different from the Xy-model describing quantum spins. Quantum fiuctuations, sur- 
viving even at T = 0, add to the thermal fiuctuations of the phases already present in the 
classical Xy-model. They tend to counteract the tendency to <Ctopological^ phase order 
which is characteristic of the 2D classical system. Taking into account only the charging 
energy with respect to ground the capacitance matrix has the form Cw = CSui and the 
influence of charging energy can be expressed in terms of the ratio 



The corresponding phase diagram (with and without dissipation) has been studied by 
various methods. Approximating the cosine phase interaction by a quadratic form in the 
phase differences corresponds to the <Cselfconsistent harmonic approximation^ [4,8,9]. One 
flnds that the critical temperature decreases with increasing a and flnally vanishes at a crit- 
ical threshold Uc on the order of 6. The same conclusion is drawn from Quantum Monte 
Carlo calculations [10]. The quantum model including dissipation is presumably also rele- 
vant for describing the superconductor-insulator transition occuring at the underdoped limit 
of cuprate superconductors [11]. It has been shown [12] that dissipation may lead to a new 
universality class characterized by an effective non-ohmic dissipation. 

The present article mainly deals with the dynamical properties of J J As. Section 2 
presents various equations of motion for the phase variables which are based on the super- 
conducting properties of the array including dissipation and charging effects. The essential 
excitations of the phase system are vortices (V) and antivortices (A) for which we derive ap- 
propriate equations of motion in section 3. The following section 4 deals with various linear 
response functions of an array, such as its electrical conductance and impedance, and relates 
such quantities to vortex variables. Various analytical and numerical theoretical techniques 




{IV) 



(1.10) 



that are currently used in order to calculate dynamic response functions are summarized in 
section 5. Effects of <Cdisorder^, ie deviations from a periodic lattice structure of the array, 
are discussed in section 6. Section 7 prepares the interpretation of experimental data by 
giving an appropriate theoretical description of the relevant measurements. The main exper- 
imental observations of the dynamic behaviour of JJAs are summarized in section 8, whereas 
section 9 presents theoretical predictions found by analytical and numerical methods. Section 
10 deals with the dynamic response of JJAs to strong applied currents. The last section gives 
some hints to other aspects of JJA that could not be covered in this article. We would like 
to warn the reader at this point that our presentation has been made to be relatively concise, 
in order to keep the article at a decent size. However, we try to give pertinent references 
where more details can easily be found. 



1.2 The Berezinski - Kosterlitz - Thouless scenario 

Treating each phase 9i as a classical variable ranging between and 2tt the Hamiltonian (1.2) 
describes the classical 2D XY-model, the phase angle giving the direction of a classical spin of 
unit length at site i. Its thermodynamic properties are governed by the Berezinski-Kosterlitz- 
Thouless [BKT) transition. The lowest energy configuration corresponds to all phases being 
equal, corresponding to the ^ferromagnetic^ ground state of the spins at T = 0. For any 
T > this long range order is replaced by <Ctopological» order due to thermal fluctuations 
[13]. Besides small amplitude variations, which give rise to <Cspin waves^ in the dynamic 
generalization of (1.2), the main finite temperature configurations are given by thermally 
excited vortices (V) and antivortices (A). Their interaction varies essentially logarithmically 
with their distance. They thus behave like a 2D Coulomb gas. In the absence of frustration 
the latter is neutral (equal number of V and A). According to the BKT scenario, V and A 
are bound in pairs below the transition temperature Tbkt- When T grows above Tbkt more 
and more pairs with large separations dissociate yielding an increasing number of <^ free ^ 
V and A. The mean distance between the latter is given by the BKT correlation length ^ 
with its particular T-dependence 

= Co exp{ J" ) (1.12) 

V-^ I -lbkt — 1 

where is a prefactor of the order of the lattice constant and h a parameter of order one. 
The topological order below Tbkt is characterized by a power-law spatial decay of phase 
(spin) correlations : 

a 

with a coupling constant J{T) renormalized by thermal fluctuations, and a finite helicity 
modulus A expressing the <Cphase stiffness^ of the system. A is given by the variation of 
the free energy of the system with respect to a ^phase twist3> : 

where the free energy 

F {a) = -keT In f d[e] e"^-^!^'"! (1.15) 



is evaluated by imposing the boundary condition 



e{la,Na) = 9{la,0) + a (1.16) 

relating the phase 9{l,m) on a given column I of the array at the lower and the upper 
boundary, respectively, the lattice sites being parametrized in units of the lattice constant 
a by {la, ma). A has a particular temperature dependence. Starting from its T = value, 
where it is given by the coupling constant J, it is a deacreasing function of temperature, since 
thermal fluctuations make the phase system softer. At Tbkt it jumps from a universal value 

Ac = -ksTBRT (1.17) 
vr 

to zero. In section 4 A will be expressed by suitable phase correlation functions and it will be 
related to the electrical conductance of the array, as well as to the vortex dielectric function. 
This allows to observe the T-dependence of A by dynamic measurements, as discussed in 
sections 7 to 9. 

Above the transition temperature the phase correlation function decays exponentially 
with the BKT correlation length given by (1.12). 

In the frustrated case there are magnetic field induced V (or A, depending on the sign of 
the field) already in the ground state. The thermal V-A pairs coexist with the latter. The 
precise nature of the corresponding phase transition in the presence of frustration is still a 
matter of debate, see ch. 8 of ref. [1]. . 

The precise link between the Xy-Hamiltonian (1.2) expressed in terms of the phases and 
the 2-d Coulomb gas can be formalized by transforming it into a vortex representation with 
the help of a dual transformation [14, 15].. This procedure allows - at least approximately 
- to separate the spin wave contribution to the partition function and yields the vortex 
Hamiltonian 

Hc = \ Y.im{R) - miR - B!){m{R!) - f) (1.18) 
RR' 

The sum runs over the plaquettes of the array, m{R) being the vorticity of plaquette R 
(in practice m = 0, 1 or -1). The interaction potential Vc depends logarithmically on the 
distance, at least for large values of the latter. The frustration / acts on the Coulomb gas of 
A and V as a <C background charge The precise critical behaviour of the BKT transition 
is obtained by introducing length scale dependent screening, which yields renormalization 
group equations [1,13,14] for the relevant parameters of the Coulomb system, namely its di- 
electric constant, related to the phase stiffness A (see section 4) and its fugacity. 

We do not go into the renormalization group approach to the BKT transition [13] since 
it concerns mainly thermodynamics. Its generalization to dynamics - our present topic - as it 
was developed by Shenoy [16] will be presented briefly in section 5.5. A very extensive review 
of the BKT transition and the Coulomb gas model has been published by Minnhagen [17]. 
Other review articles, such as [18,19], are more directly concerned with dynamic properties 
of J, J As. In order to interpret correctly measurements on J, J As in the light of the BKT 
scenario one should, however, bear in mind that the latter (as critical phenomena in general) 
is in principle valid for the thermodynamic limit of an infinitely large system. In section 



7.1 we will mention finite size effects that can alter the picture to some extent, in particular 
through the presence of unbound vortices below the transition temperature. 



2 The Dynamics of JJA 

2.1 Equation of motion for the phases 

In order to derive equations of motion for a JJA we start from our Hamiltonian (1.9) involving 
Josephson coupling and electrostatic energy. It is useful for later purposes to include a possible 
coupling of the array to external currents Ii flowing into (or out of) site I 

Hi = J2lA (2-1) 

The resulting canonical equations read 

0i = ^ = {2efT.(C-')ii'Pi' (2.2) 
on ^, 

^' = -1^ = - E -^f' «"^(^' - ^'0 - Ii (2.3) 

Due to (1.6) and (1.8) the second equation is the continuity relation expressing charge 
conservation at each site of the array, the right hand side being the sum of supercurrents and 
external current flowing into site I. 



Dissipation, due to normal current between islands and from islands to ground, is now 
easily added to the charge conservation expressed by (2.3) which then reads : 

Pi = -Y. Jn' s\n{ei - Ov) - J2iR-')a'i^f{0i - Ov) - ^oH^fOi - h (2.4) 

Here we have introduced the matrix Ruf of junction resistances and the resistance Rq to 
ground, and we have used the Josephson relation between the time derivative of the phase and 
a corresponding voltage. Equations (2.2) and (2.4) correspond to a description of the array 
in terms of an effective circuit involving capacitive, resistive and Josephson (^ inductive ^) 
elements. Neglecting the <C local 3> charging energy and normal current (ie between a given 
site and the ground) yields the so-called <C Resistively and capacitively shunted junction 
{RCSJ) model >. 

Capacitive effects are important for ultrasmall superconducting islands. The interplay 
between charge and phase fluctuations lowers the superconducting transition temperature 
and can drive it completely to zero, leading thus to a quantum superconductor-insulator 
transition, as discussed in section 1.1. We will mostly concentrate on dynamic phenomena 
of <^ classical arrays ^ for which charging effects are unimportant. This corresponds to a 
motion without inertial acceleration and thus equation (2.4) with the left hand equal to zero 
will determine our phase dynamics. The two most frequently used equations of motion are 
then distinguished according to the type of phase damping involved, corresponding to the 
type of normal current dissipation : 



a. Resistively shunted junction - RSJ - dynamics 

Only the ^ bond damping ^ (the second term in (2.4)) is kept, describing dissipation by 
currents flowing through the junctions. This dynamics is therefore also called <C total current 
conserved With Rj being the junction resistance the RSJ equation for the phases then 
reads : 

A-jrp 1 

V I" 

Here we have introduced the Green function of the lattice which allows to write the sum over 
bond differences of the phase as a sum over sites : 

V I' 

b. Time dependent Landau-Ginzburg - TDGL - dynamics 

Taking into account only site damping due to losses to the ground (only Rq is considered 
in (2.4)) yields the simpler <Crelaxational3> equation of motion 




This corresponds to the usual time dependent Landau-Ginzburg equation representing <^ 
nonconserved 2> dynamics (charges are ^ lost 3> by currents going to the ground). 

In the presence of an electromagnetic field the phase differences occurring in the various 
current expressions have to be replaced by the corresponding <C gauge invariant ^ combi- 
nations with the applied vector potential, as indicated in (1.4). 

Dynamic equations, such as (2.6) or (2.7) can now be used for various analytical consid- 
erations, such as the derivation of an equation of motion for the center of a vortex (section 
3) or the calculation of dynamic linear response functions (section 4). They are also at the 
basis of numerical simulations of equilibrium dynamics (section 5.4) and of the response of 
the array to very strong driving currents (section 10). 

2.2 Microscopic derivation of the coupling constants 

The various coupling constants (resistances, capacities and Josephson couplings) of the RCSJ 
model describe the physical properties of single Josephson junctions. They can be derived 
from a microscopic transfer Hamiltonian, describing the coupling between the two super- 
conducting islands through a normal metal or an insulator by a tunneling interaction and 
including Coulomb repulsion between electrons. Besides numerous earlier publications on the 
physics of Josephson junctions useful reviews, focussing in particular on very small junctions 
which can now be fabricated, have been given by Schon and Zaikin [20] and by Simanek 
[4]. Using either an imaginary time or a real time representation the quantum mechanical 
partition function can be rewritten in terms of an action involving the phase difference be- 
tween the two superconductors. Its dynamics is indeed governed by a capacitance and by 
generalized dissipation and Josephson coupling terms. 



(2.6) 



3 Vortex dynamics 



3.1 Classical equations of motion 

Vortex excitations are the main elements determining the physical properties of a J J A, in 
particular in connection with the Berezinski-Kosterlitz-Thouless scenario, as explained in 
section 1.2. Thus it is important to account for them explicitly in dynamic calculations. A 
classical equation of motion for a vortex configuration can be found, in analogy with the way 
the dynamics of localized excitations in magnetic models is treated [21]. We start from one of 
our general equations of motion for the phases of the array derived in section 2.1 which can 
be transformed into an equation for the bond supercurrent given by (1.6). In order to prepare 
a continuum description of the array it is useful to number each site (of a square array) by 
an integer / and to attach a horizontal and a vertical bond to each site, distinguished by the 
index s = x or y, respectively(see figure 1). The bond current Kn^ (1-6) is then replaced by 
Ks{I) and the phase difference (1.4) by (ps{I)- 



Figure 1 : Parametrization of a square lattice : each site is labelled by an integer /. 
Two nearest neighbor bonds attached to each site yield the two cartesian directions x and y, 
denotes by the corresponding unit vectors pointing along the plaquette /. They determine 
the X and y components of vector and tensor quantities in equations (3.1) to (3.5). 

Our equations of motion for the phases, derived in section 2, can be used for determining 
the time evolution of Ks 

^"—Q^ = T^U^sAi - i')KAi') + vsAi - i')Mi')} (3.1) 

s'l' 

where Es{I) = cos(j)s{I) and the ^dynamical matrix^ M ensures the correct coupling be- 
tween bonds. 

In (3.1) we have specialized to the case of a local capacity, but we still admit both types 
of damping discussed in section 2.1 : 

, _ J Sss'Sii'VL for local damping , , 

"^'"'^ \ Mss'il - I')VT for bond damping ^'^ ' 

The friction constants are proportional to the inverse resistances showing up in (2.4), ie the 
ground (local damping), respectively the bond resistance (bond damping). 



A moving vortex is now introduced by a trial solution [22] 



Ml,t) = cp'f\l,t) + ips{I,t) (3.3) 

where 



{I,t) = ar c sin Ki''\l\K{t)) (3.4) 



The first term represents the circular symmetric <Chedge hog^ phase pattern with a time 
dependent center R(t) and the second contribution describes deviations from this form. The 
form (3.2) is substituted into the phase equation of motion (3.1) yielding an equation for the 
bond currents K^"^^ forming the vortex pattern 

Co ^/ ' =-E'{i)Y,{JM.,Ai-i')[EAi>Ai')+K^':\i')]+VsA 

s'l' 

(3.5) 

(v) 

We have neglected the small deformation field (p wherever it multiplies terms involving (f>s . 
The first term on the right hand side of (3.5), involving the dynamical matrix M, represents 
the Peierls force acting on the moving vortex which originates in the discrete lattice structure 
of the J J A. Indeed, the potential energy of the phase pattern of the vortex increases when 
its geometrical center moves away from the middle of a plaquette. It reaches its maximum 
value when the vortex center crosses the bond separating two adjacent plaquettes. The height 
Eb of this energy barrier depends on the lattice structure : Eb « 0.199 J for the square and 
Eb ~ 0.043 J for the triangular array [23]. Going over from the discrete lattice to a continuum 
and integrating out the phase variables in a suitable way (see [21] and [22] for more details) 
allows to obtain a Newtonian equation for the center R of the vortex 

M^il = -TR - 2TiqaJz x (p (3.6) 

The parameter q = ±1 determines the sign of the excitation (V or A). The vortex mass My 
is expressed in terms of the grain capacity Co and a <^ core radius ^ 6 of the vortex patterns 
(which can be taken of the order ot the lattice constant a) : 

M. = Co'^i-? (3.7) 

Obviously, in the absence of an underlying discrete lattice structure there is no more energy 
barrier for the moving vortex. The two types (3.2) of dissipation, although yielding different 
equations for the phase field, lead to the same friction term in (3.6), the respective friction 
constants being only quantitatively different : 



VL^^i^)^ for local damping 
^^T"(a)^ for bond damping 



(3.8) 



One should stress here that different calculations may yield different results for friction con- 
stant and mass: <Csimple^ calculations like the above skteched trial function approach yield 
a <SbareS> mass, whereas taking into account the coupling of the vortex configuration with 
thermal fluctuations introduce a renormalization to mass and friction, which can even be 
(logarithmically) divergent due to the 2D ^substrate3>, see, for example, [24]. The second 
term on the right hand side of (3.6) can be interpreted in different ways : it can represent 
the force acting on the vortex when a supercurrent, given by a phase gradient ip = V9, is 
flowing through the sample. It is usually called "Lorentz force" [1], [25], in analogy to the 



force a current-carrying wire experiences in a perpendicular magnetic field. In our later cal- 
culations we will insert into equation (3.5) the phase gradient VO resulting from the presence 
from another vortex located at Ri in order to introduce the effect of V-V interaction in the 
equation of motion : 

M^R = -m + gigsflV^ ~^\. (3.9) 

[ti. - ±tij^ 

When capacitive effects in the array can be neglected (for sufficiently large capacities or for 
low enough frequencies) the mass term disappears from (3.6). 



3.2 Quantum dynamics 

At sufficiently low temperatures the JJA, characterized by the Hamiltonian (1.9) and an 
appropriate dissipation function, should be treated by quantum mechanics. As sketched in 
section 1.1 this is necessary when the <Cmass term^ in (1.9) becomes comparable to the 
Josephson coupling, ie when the <Cquanticity ratio^> a given by (1.11) is large enough. 

The corresponding dynamics of quantum vortices has been studied on the basis of a 
quantum action for the phases, which allows to treat dissipation on the same footing as the 
Hamiltonian terms. Vortex variables are introduced in two ways. Generalizing the Villain 
transformation, used for the same purpose for the classical XK-model [14] leads to dynamic 
plaquette vorticities [26]. The transformation is performed on the quantum form of the par- 
tition function involving an action including charging energy and dissipation. As a result one 
obtains a new action in which the (true) electric charges interact with V and A, acting like 
^topological charges^. Integrating out the electrical charges yields an effective Lagrangian 
for the V-A system that attributes effective mass and friction to the latter. Alternatively, the 
action is expressed in terms of vortex positions and momenta which are again introduced by 
a trial form for the dynamic phase field of a vortex [27], as it was done above for the classical 
case. 

Quantum vortices having a non-zero mass show interesting quantum effect [25], in partic- 
ular when the lattice pinning potential is taken into account. While classical particles would 
oscillate in the minima of this potential and jump across barriers by thermal excitation, quan- 
tum vortices can tunnel through pinning barriers. Moreover, the quantum states of a vortex 
in the lattice potential will be a Bloch state with some given momentum. As in the case of 
band electrons Bloch oscillations can occur, as the particles move across the Brillouin zone, 
being driven by en effective electric field. Quantum interference and quantum localization in 
a random pinning potential (see also at the end of subsection 5.3) are another manifestation 
of the quantum nature of vortices. 

We do not give more details here, since our subsequent analysis will mainly be focussed 
on <C classical ^ arrays. 



4 Linear response of a JJA 



4.1 Kubo formalism for conductance and impedance 



In this section we prepare the theoretical description of experiments that ^ test ^ the dy- 
namic behaviour of a J J A. This can be done in different ways by applying suitable external 
disturbances to the system. 

The frequency dependent conductance G{uj) is obtained from measuring the current that 
flows in the array as a response to a small time dependent external voltage. Alternatively, 
when a weak external current is fed into the array the impedance Z{lo) determines the result- 
ing voltage drop across it. Explicit expressions for G{io) and Z{uj) in terms of suitable time 
dependent correlation functions are obtained in this subsection in the framework of linear 
response theroy. In subsection 4.2 the response functions G{lo) and Z{uj) will be expressed in 
terms of quantities characterizing the vortex system. Other experimental approaches, such 
as the two-coil technique or the measurement of flux noise, will be described in section 7. 

Since the JJA dynamics involves dissipation it is convenient to interpret the equations 
of motion (2.3, 2.5) as Langevin equations for the phase and the conjugate momentum by 
adding a random force. For calculating response functions it is useful to go over to the 
corresponding Fokker-Planck (FP) equation [38,39], governing the probability distribution 
p{6i, -Pi, ..) of the phases and the momenta : 

^ = LfpP (4.1) 

The FP-operator Lpp = L + L^x has a contribution L containing convective and dissipative 
terms present for an isolated array and Lg^ describing the coupling to external sources. 
Besides the coupling to external currents /; in (2.1) an external electric field can be applied 
through a time dependent vector potential Ani in the phase differences (1.4) yielding the 
following form of L^x 



t 1 1 t £ 



Each term in (4.2) has the general structure L^x = X * F{t) where X is an observable of the 
array expressed in terms of the angles 9i and the conjugate momenta Pi and F{t) is a time 
dependent external force F{t). In the FP formalism the linear response of an observable O in 
the presence of such an external perturbation is given by splitting the probability distribution 

P = Po + Pi (4.3) 

into the canonical equilibrium part po and the deviation pi due to F. They are, respectively, 
given by 

Po = ^exp{-(3H) 



Pi{t)=t^dTe^'--)^L,^{T)po 



(4.4) 



assuming that the perturbation is adiabatically switched on at t = — oo. The expectation 
value of O is then given by 

{0{t)) = f dr R{t - t)F{t) + Oo (4.5) 

J —oo 



where Oo is the equihbrium expectation value of O. The frequency Fourier transform of 
the response function R is given by the retarded time dependent correlation function of the 
variables O and X. 

poo 

R{lo) = / dt e'^^^TriOe^^Xpo) (4.6) 
Jo 

the symbol ^ Tr ^ meaning integration over phase and momentum variables. 

We now apply this general scheme to the two experimental configurations mentioned 
above : 



a) For calculating the conductance G{uj) the observable O is the current (supercurrent 
plus normal current) flowing through some junction 

O = Jiv s\Tv(t>w + {R^^)w{^fH' (4.7) 

47re 

The action of L^x is given by the second and the third term of (4.2) which, according to (4.4), 
will act on po- The derivative with respect to Pi is easily evaluated using (4.4), (1.9) and 
(2.2) : 

^ = -fjp,Y^{C'%.Pv = -UpoQi (4.8) 
dPi ^ 

We specialize to the situation where the applied electric field is homogeneously spread 
across the Nr^ sites of the array, say in rr-direction. Thus the time derivative ot the vector 
potential applied to a bond going from site I to site I' in the (positive) x-direction is given by 

Aw = = -Am (4.9) 

As for calculating the static helicity modulus it is important to impose suitable boundary 
conditions on the phases for the conductance calculation. Since we deal with the variables 
Q\ and Pi in the present context we use periodic boundary conditions {PBC) for the phases. 
Later on, in section 5.4 on numerical simulations, boundary conditions will be introduced 
which are more appropriate when the relevant variables are vortex positions. Taking (4.9) 
and the PBC into account sums like 

s=Y.{ei- Ov) (4.10) 

extending over all bonds in x-direction vanish. Thus the total normal current in (4.7) reduces 
to 

4"^ = -^E^«' = ^^- (4-11) 

(«'> 

which gives the ohmic contribution to the conductance. Moreover, for a regular array the 
parameters Jw and rjui are equal for all bonds (//') so that, due again to (4.9) and the PBC 
the third term in (4.2), representing the coupling of the applied vector potential in the normal 
current, also vanishes. We only have to evaluate the response of the supercurrent in (4.7) to 
the second term of the external perturbation in (4.2). It has two contributions : 

J2 rt 

{K^{t)) = J ^{cos(t)w)Aw{t) - 13— j (sin <?!'//'(*) sin (r)(/)„„'(r))^„„/(r) 

(4.12) 



Partial integration of the second term and Fourier transforming with respect to time leads 
to the final result 

{K,{lo)) = Gs{oj)E^ (4.13) 
where the superconducting part of the array conductance 

Gs{oj) = — - PF,{uj) (4.14) 

is expressed in terms of the helicity modulus 

J2 

A = J ^ {cos (pw) - (3— ^ (sin 0;/' sin (4.15) 

ill'):, '-^ {lV)x {nn')^ 

and a dynamic current-current correlation function 

F,{oj)=j dte'^'j^Y.T.^''''''t>iv{t)sin(t>nn'm (4.16) 

^ {ll')x {nn')x 

Finally, the total conductance G{uj) is given by 

G{Lo) = ^-(3Fe{u;)+7] (4.17) 

The superconducting contribution to G has already been derived by Shenoy [28]. 

Thus the dynamic conductance G is expressed in terms of an effective circuit with parallel 
superconducting and normal channels. This picture can be made more explicit by writing G 
as 

iloL[lo) K\uj) 

with two real, frequency dependent response functions : an effective inductance L and an ef- 
fective resistance R. Comparison of (4.17) with (4.18) shows that the helicity modulus which 
is one of the basic features of the BKT transition can be extracted from the low frequency 
behaviour of the dynamic conductance. 

b) The impedance Z{ijo) is obtained by calculating the expectation value of the voltage 

Viv = ^{9i-9v) (4.19) 

across a given junction under the influence of an applied current acting through the first 
term of L^x in (4.2). In order to allow the external current to enter the array open boundary 
conditions have to be used for this calculation. Again using (4.8) for the derivative of po in 
(4.3), the response function is immediately found to be 

{Vu'{t))=f dTY,mt)-9,{t))9nir))In{T) (4.20) 

The full impedance is found by summing the voltage drops over all bonds across the 
sample in x-direction and letting an equal current / flowing in on the left side through the 
sites (0,ma), respectively out on the right sides {Nxa,ma) of the array. This yields the final 
expression for the impedance 



^(^) = TT E / dt e'^\{9i{t) - OvitWo^ma - 9n„)) (4.21) 



The expressions for G and Z in terms of the corresponding dynamic correlation functions 
do not seem to have much in common. It is thus legitimate to ask the question whether the 
expected relation 

G{u:) = Z{uj)-^ (4.22) 

holds. Expressing these correlation functions in the following subsection, for a classical array, 
by suitable vortex variables will indeed confirm (4.22). Alternatively, expressing the partition 
function of an array including capacitive effects and damping terms as an imaginary time 
functional integral, and going to a dual representation using the electrical charge on each site 
rather than the phase as variables, Panykov and Zaikin [29] have shown that the correspond- 
ing linear response expressions for G and Z indeed satisfy (4.22). However, relation (4.22), 
as plausible it may look, need not necessarily hold. Indeed, it has been shown [29,30] that 
for a quantum array the linear voltage response to an external current is not equivalent to 
the inverse of the linear current response to an applied voltage. The difference between the 
two arises because the voltage bias suppresses some important quantum fluctuations of the 
phase (instantons) due to the fact that the total phase difference across the sample is fixed by 
the applied voltage drop. Therefore the phases are hindered in their quantum fluctuations. 
The latter are not suppressed by an applied current since open boundary conditions are then 
used. As a result, the impedance Z sees the sample <C less superconducting ^ than the 
conductance ! The same effect could, in principle, also occur with thermal fluctuations in the 
classical regime, but an explicit example seems not to have been constructed. 



4.2 Dynamic correlation functions and dielectric function for the vortex 
system 

As stated in section 1.2 our vortex-antivortex system is modelled by the 2D neutral Coulomb 
gas (2DCG). It is therefore useful to introduce some concepts that are currently used when 
treating the dynamics of the 2DCG (see, for example, the chapter ^Charged fluids^ of ref. 
[31]). For a hydrodynamic description, adapted to treat various <Cmacroscopic» experiments 
treated in section 7, the basic variables are charge density p{r) and current density j(r), the 
spatial Fourier transforms of which are related to the positions R; and the sign = ±1 of 
the charge of the N particles by 

(4.23) 

In section 5.3 we will use Mori's formalism in order to calculate dynamic correlation 
functions of these basic observables, such as the <C charge density correlator^ 

roo 

p,z)= / dt e''\p{p,t)p* {p,0)) (4.24) 
Jo 

where the average is over a canonical ensemble and the time dependent particle positions 
have to be obtained by integrating suitable equations of motion. The equal time correlation 
function 

Scip) = {\p{p)\') (4.25) 
is the charge structure factor of the 2DCG. (j)pp and Sc are related to the dynamic charge 
susceptibility x by [31] 

X(P, z) = --i-[5,(p) - z^ppip, z)] (4.26) 



which, in turn, is hnked with the dynamic dielectric function e{uj) by 



27r qln 



^ = lnn[l + ^x(p,-^-)] (4.27) 

where Qo is the charge of the particles and n their density. As usual in electrodynamics the 
dielectric function can be expressed in terms of a dynamic conductance o"„(a;): 

e{u;) = 1 - 27T^^ (4.28) 

Here we have used the letter a and a subscript v in order to distinguish this ^ vortex 
conductance related to the mobility of the vortices and thus to their response to an ap- 
plied force, from the physical electric conductance G of the J J A. 



For many purposes - in particular in order to illustrate the relation G{u>) = Z{u>)^^- it 
is now useful to express our response functions directly by these dynamic vortex correlation 
functions. Various procedures have been proposed in order to transform the XY problem 
from its phase form (expression (1.2)) into a vortex representation. An approximate dual 
transformation allows to separate linear (spin wave) and vortex excitations, yielding the 
vortex Hamiltonian (1.18). This is particularly useful for dealing with equilibrium thermo- 
dynamics, but the same approach has also been applied to the action in a functional integral 
representation of the quantum version of the XY-model. 



Another approach yielding insight in the dynamics of the system [33] consists in approx- 
imating the lattice phase Hamiltonian by a an energy functional for continuous phase and 
momentum fields 

H = ij d\{Vef + dVP(r) (4.29) 

and separating the phase field 6 into a <ti regular S> (spin wave) part (/) and a vortex part ■0, 
given by 

61 = + ^; V X V0 = 0, V • V0 = (4.30) 

Cote and Griffin [33] transform the equations of motion for the phase and the conjugate 
momentum fields, without dissipative terms, into Maxwell equations in which the vortices 
act like effective charges being the source of effective electric and magnetic fields related to 
the phase field by 

E = V2^ve X z 

(4.31) 

z being a unit vector perpendicular to the array. The charge of a vortex is related to the 
original coupling constant in (4.29) by 



ql = 27rJ (4.32) 

This concept can easily be generalized to include dissipation by adding a local damping term 
(normal currents flowing to the ground) to the equation of motion as in section 2. Faraday's 
equation is then generalized to 



(4.33) 



with 

(4.34) 

being the <^ speed of hght ^ in this effective electrodynamics of the array. Expressing super 
and normal current densities by 

K, = jve 

(4.35) 

allows now to express the response function found above for G and Z in terms of vortex 
quantities. 

For the impedance Z we have to establish a linear relation between the total current, 
I = + In, and the voltage across the sample. We start from Maxwell's equation 

1 ■ 2-7T 

V X B = -E j„ (4.36) 

c c 

being the vortex current density and rewrite it in terms of the time derivative of the phase 
gradient 

V^ = ^=i=(lxE-27rixj^) (4.37) 

Using (4.31) and (4.35) we replace the effective electric field by Kg = I — Kn- Then we relate, 
within our <C effective electrodynamics^, the vortex current density to the driving electric 
field through the vortex conductance (T„ : 

j^, = tr^E (4.38) 

which yields 

{l + zuj-j -r]—^)V9 = I{— —) (4.39) 

The dynamic vortex dielectric function is introduced by (4.28). Using the dynamic Joseph- 
son relation 

V = —V9 (4.40) 

2e 

we finally express the impedance Z = y in terms of ey{uj): 

Z(.) = |— ^ (4.41) 

For the dynamic conductance, given by (4.17), we first note that the helicity modulus A 
is directly related to the static vortex dielectric function[15]: 

A = -r^, (4-42) 

ev[uj = 0) 



In the continuum description the supercurrent density is given by (4.35). For our PBC 
the regular part of the phase gradient does not contribute to the dynamic current-current 
correlator in (4.16). According to (4.31) the singular part stands for the effective electric field 



related to the vortex charge density by Poisson's equation. Using also = 2ttJ this allows 
to relate (4.16) to the dynamic vortex charge correlator (4.24) : 



Fe(a;) = pj\2^fn. lim '^^^^P', (4.43) 

riv being the vortex number density. The way towards the vortex dielectric function is found 
by using the relations (4.26) and (4.27), once for finite frequency oo and once for = 0, which 
immediately yields 

F,(a;) = -[^-4iT] (4-44) 
luJ e„(w) e„(0) 

Substituting this result and (4.17) in the expression (4.20) for G leads to 

G(a;) = — ^ + r? (4.45) 

Manifestly it is the reciprocal of the impedance given by (4.41). Expression (4.45) is a partic- 
ular case of the form (4.18), in which the original bond resistance becomes frequency depen- 
dent through the vortex dielectric function. Expressing the latter by the vortex conductivity 
through (4.45) one obtains 

Giu) = + 77 (4.46) 

Thus, the superconducting contribution to G is directly related to the vortex mobility : 
G(0) is inversely proportional to the latter. In particular, when the vortex system has a non 
vanishing mobility G(0) is finite and superconductivity is destroyed. 

5 Theoretical techniques for calculating dynamic vortex re- 
sponse function 

5.1 Early approaches 

Ambegaokar et al[34] calculated the dielectric function e{uj) of the vortex system 

e{uo) = l + 2^x{^) (5.1) 

for independent bound pairs by averaging the dynamic susceptibility x of a dipole of length 
r over a Boltzmann probability distribution y3(r, T) = e"^^^"^'''^ of finding a V-A-pair of 
length r at a given temperature. A simple form for the effective interaction V (the Coulomb 
interaction devided by a dielectric constant) yields 

^(r) oc (r)-4^ (5.2) 
a 

Alternatively one can use a length dependent static dielectric function [34 and 35]. The 
resulting dielectric function 

e{uj) = 1 + ^ r rdr ^^^^-r (5-3) 

has a power law dependence on frequency with exponents that depend on temperature. 



For the free vortex system, above Tbkt, the simplest form for e is given by Drude's form 
obtained by inserting a frequency independent vortex conductivity cTo which is proportional to 
the inverse of the friction function in the equation of motion (3.6). This yields the well-known 
simple form for the inverse of the dielectric function : 

Re(^) = 

(5.4) 



5.2 Minnhagen's phenomenology (MP) 

This approach [17] starts from a Debye form of the wave number dependent static dielectric 
function = p^pfc? with a screening wave vector kg and goes to Re{l/e{uj)) by replacing 
the wave number k by frequency uj through 

e = ^ (5.5) 

which corresponds to taking into account diffusive motion of the particles with diffusivity D. 
Performing a Kramers-Kronig transformation one obtains 

(5.6) 

Since e{uj) = 1 + a{uj)/{—iuj) this amounts to using a vortex conductivity with a peculiar 
non-analytic frequency dependence for small uj 

i?ea(a;) = -^(lnM)-i 

(5.7) 

Ima{u:) = ^(In^) 

At first sight one would tend to reject this result, since - had one made the more correct 
replacement 

u2 _ 

as one should for a diffusive dynamics - one would come back to Drude's form rather than to 
(5.6) and (5.7). It turns out, however, that simulations and experimental data confirm the 
MP form of the dynamic response in an appropriate frequency domain, as we will show in 
sections 8 and 9. 



A simple criterion allows to distinguish the MP form (5.6) from the more simple Drude 
expression : the ^peak ration 

/m(^) 

between imaginary and real parts of the inverse dielectric function, evaluated at the frequency 
ujo where ImiXje) has its maximum, is given by r = 1 for Drude's form and by r = 2/7r for 
MP [17]. 



5.3 Equations of motion (Mori, Fokker-Planck) 

The dynamic charge density correlator (f)pp(p,z) defined in (4.24) is a central quantity for 
treating vortex dynamics. Indeed, it is directly related to experimentally accessible quantities 
such as the dynamic dielectric function (see section 4.2) and the flux noise spectrum that 
will be introduced in section 7.3. It is therefore important to rely on systematic methods for 
calculating such classical time dependent correlation functions. 

Mori's approach [36] for evaluating dynamic correlation functions of systems for which 
relevant static correlation functions are known is particularly useful for our 2D Coulomb gas. 
It can be based on the two hydrodynamic variables p and j introduced in section 4.2. The 
time evolution of these variables is easily found by using the equations of motion for the 
center of a vortex (see section 3) including force and friction terms. Mori's technique yields 
a formally exact expression for (ppp{p,z) 

Mp,^) = —^^^^ (5-9) 

where Sc{p) is the charge structure factor (4.25) of the vortex system and 7(p, z) - the ^ 
memory function ^ given by the dynamic correlator of the ^random forced - summarizes 
the complicated effects of the Coulomb interaction on charge dynamics [35] . Various approx- 
imations can be used in order to obtain manageable expressions for 'y(p,z). We first focus 
on the region above the BKT temperature treating the vortex system in its <Cmetallic^ 
phase, ie by neglecting possible binding in pairs. We assume that the zero wave vector limit 
7(0, z) = j{z) is sufficient to describe hydrodynamic phenomena. The general form of this 
quantity is given by [35] : 

7(^)^r+ ^° 3 ^k-kVe(k)y;(kO / dte-^*(n*(k,t)p(k,t)n(k',0)Kk',0)) (5.10) 
2fcBTnt,iF ^ Jo 

Here Q is the total area of the system, Vc{p) = is the Fourier transform of the 

Coulomb potential and n(p) is the particle density 

n(p) = ^Ve^P-i^' 

corresponding to the charge density in (4.8). In the above statement we have - as it is usually 
done - neglected the action of the <CMori projector^ [36] in the time evolution of the higher 
order correlation function on the right hand side of (5.10). In the absence of interparticle 
interaction 7(2) simply reduces to the bare friction T. Thus j{z) is a generalized friction 
function and its inverse is proportional to the <C charge mobility^ of the system. This can 
be seen by inserting the specific form (5.8) of (j)pp{p,z) into the relations (4.26) and (4.27) 
determining the dielectric function. One finds the form (4.28) with 

2Trqlny 

^v{uj) = , . . (5.11) 
7(-zu;) 



This expressions illustrates the role of the inverse of the friction function as a ^charge 
mobility^. 



As usual in treating interacting systems a higher order correlator, involving both charge 
and number densities, shows up in (5.10). In the spirit of the <C mode coupling ^ approxi- 
mation, currently used in liquid dynamics [37], this higher order correlator is factorized into 
a product of (j)pp{p, z) and of the corresponding number density correlator 

POO 

(l)nn{p,z)= dt e-''{n{pt)n*{pO)) (5.12) 
Jo 

However, in order that such a factorization be a useful approximation, the long range 
Coulomb potential - which induces long range correlations between particles - should be 
replaced by the screened interaction, as it is customary in treating the dynamic properties of 
Coulomb systems [40,41] 

"""^"^ = Mi} 

which, in this simple form, involves the Debye screening wave number 

P} = (5.14) 

determined by the charge Qo and the density Uf of free vortices. For the static structure 
factor Sc{p) we use the simple form 

2 

Scip) = (5-15) 

V +Po 

with 

2 kBT 



Po o 9 

ZTTq^n. 



V 



which satisfies the sum rule [31] for the metallic phase of a Coulomb gas. Note, that po, the 
relevant wave number for Sc{p), involves the total particle density. For the number density 
correlator (5.11) we use a simple diffusion form 

Tl 

(l)nn{P,z) = (5.16) 



Z + 



Equation (5.9) is still an implicit expression for 7(2) since the right hand side involves 
the same quantity again through (j)pp{p,z). To lowest order in the interaction 7(2) can be 
evaluated by replacing it by F under the integral on the right hand side of (5.9). Results for 
7(2;) and thus for the flux noise spectrum and the dielectric function, obtained by using this 
procedure will be mentioned in section 9. 



Mori's approach also allows to take into account pinning (Peierls) forces acting on the 
vortices, as they result from equation (3.5) (see the discussion following that equation). One 
then obtains another contribution to the friction function 7(2:), having a similar structure 
as (5.10), with the Coulomb force being replaced by the pinning force. This approach has 
been applied [42] to take into account the effect of random pinning in disordered lattices 
(described in section 6.). A mean-field like description of disorder induced localization of 
vortices is found in this way [42]. 



The Fokker-Planck equation, presented in section 4.1 in connection with hnear response 
calculation using phases as variables, can also be used for determining the vortex distribu- 
tion function, providing thus another systematic approach for evaluating dynamic correlation 
functions. It has been used by Timm [43] for calculating the flux noise spectrum in an inde- 
pendent pair approximation. The corresponding Fokker-Planck operator takes into account 
the free diffusion of pairs, the internal interaction between the partners and a possible anni- 
hilation of the V and A constituting the pair. The results of this approach will be mentioned 
in section 9. 

5.4 Numerical simulations 

The dynamic properties of JJAs have been studied by solving numerically different types of 
equations of motion: 

a) for the phases of the XY-model using RSJ dynamics (see our equation (2.5)) : Refer- 
ences 44-54. 

b) for the phases of the XY-model using TDGL dynamics (see our equation (2.7)) : Ref- 
erences 44-46, 50, 53, 55-60. 

c) for the positions of the charges in the 2D V-A-Coulomg gas (see our equation (3.9)) : 
References 59, 61-63. 

d) for the 2D lattice Coulomb gas using Monte Carlo dynamics : References 64, 65. 

Details about the simulation technique can be found in these references. The correct 
choice of boundary conditions is particularly important. Indeed, various quantities that are 
relevant for the BKT transition, such as the phase stiffness, describe the response of the 
system to forces applied to its boundaries. For simulations of the Xy-model (a and b) 
periodic boundary conditions (PBC) imposed in the phase angle 9i are most commonly used 
when the total current in the sample is zero (see our discussion of boundary conditions in 
the context of linear response calculations in section 4.1). However, as shown by Olsson [66], 
PBC for the phase angles lead to non-periodic boundary conditions for the vortex interaction. 
On the other hand, Coulomb gas simulations (c and d) prefer PBC for the motion of the 
point particles representing the vortices. This requires imposing <Cfluctuating twist boundary 
conditions^ (FTBC) on the phase angles [66]. In this case a variable D = {D^, Dy) is added 
to the phase difference (1.4) between neighboring sites / and /' : 

4^11' ^ 4^11' - ■'D (5.17) 

the vector D being projected onto the unit vector e of the bond {II')- For the phase angles 
themselves PBC are used. D(t) is itself a dynamic quantity which is related to vortex motion 
: whenever a vortex has crossed the sample of length L from one side to the other D changes 
by 27r/L. Its equation of motion is established by considering [45] current conservation across 
the sample 

^ + J ^ s\n{4>w -D^)+ij^ = (5.18) 

Here, the first term expresses the normal current in x-direction in terms of the voltage 
drop Vx and the normal resistance R, and the last term represents the thermal noise current 



used in the usual Langevin equation in order to impose a given temperature T on the system. 
By the dynamic Josephson relation the voltage is proportional to the time derivative of the 
total phase change across the sample. Since PBC are imposed on 9i it is directly related to 
the time derivative of ■ 

V, = -^LD, (5.19) 

Combining these relations leads to an equation of motion for D [66] . By considering the 
Fokker-Planck equation for the probability distribution of the phase angles and the additional 
variables {Dx,Dy) one obtains the standard white noise correlations for the noise current : 

2T 

{Vx{t)Vxm = jjS{t) (5.20) 

A current passing through the sample can either be incorporated in the current conser- 
vation condition that determines D(t), or - alternatively - open boundary conditions for the 
phase angles can be used (see section 4.1). The former approach has the advantage that one 
can continue to use periodic boundary conditions, which is more convenient for numerical 
calculations. 

5.5 Dynamic scaling 

The BKT scenario governs the phase transition in a 2D model involving a complex order 
parameter with an amplitude and a phase. It has the particularity that the correlation length, 
diverging as usual when T goes to its critical value from above, stays infinite for all T < Tkt- 
Thus the usual dynamic scaling relation 

F{k,Lo) = F{k^,ujn 

(5.21) 

F{r,t) = F{r/C,tm 

for functions depending on space and time variables, is valid a priori only above Tbkt- The 
dynamic critical exponent z summarizes the essential information about vortex dynamics. 
Numerous calculations, mainly numerical, have been aimed at obtaining z. Apriori, one ex- 
pects z = 2 [50], according to the usual classification of dynamic critical phenomena [67], 
which is indeed bourne out by various calculations, and it is also reflected by the flux noise 
measurements of ref [68]. It corresponds to the fact that motion of vortices is of diffusive 
nature. 

Below Tbkt, where the model remains critical, the situation is more complex. Since 
^ = oo the relevant scaling now involves other typical lengths, such as the sample size or 
the size of the detection device (such as the coil used for flux noise measurement, see section 
7). Morevoer, z depends on the dynamic phenomenon to which it is applied and thus on the 
boundray conditions used. Jenssen et al [50] give a rather complete overview over previous 
results and over their own calculations. They have used RSJ, as well as LG dynamics for the 
phases (see section 2), applying, for both cases, either periodic or fluctuating twist boundary 
conditions and they have analysed resistance, supercurrent scaling and the short time relax- 
ation of a non-equilibrium configuration. 

An interesting observation has recently been made for the three dimensional(3D) XY- 
model : there is evidence [59] for the existence of two dynamic critical exponents signalling 



the dependence of dynamic correlation functions on wave number k and on the correlation 
length ^: 

F{k, t) = F{kC, tC"" , tk-') (5.22) 

with Zo = 1.5 and z = 1, respectively z = 2 for the 3D lattice Coulomb gas, respectively the 
3D Xy-model. 

Prom an experimental point of view a particularly interesting application relates the value 
of z below Tbkt to the non-linear current- voltage characteristics V ~ (see section 7.1). 
Ambegaokar et al [34] found 

a{T) = ^ + 1 (5-23) 

where J{T) = J/e{T) is the effective XY-coupling constant, renormalized by spin wave and 
VA-pair fluctuations represented by the static dielectric constant e(r). The non-linear I- 
V relation relies on the fact that VA pairs that are separated by an applied current (and 
thus contribute to a finite voltage like ^ free ^ vortices) can recombine, see section 7.1. 
Expression (5.22) is found be setting the pair breaking and the recombination rates equal to 
each other [34]. At the critical temperature one obtains a = 3, which together with z = 2 
satisfies the dynamical scaling relation [69] a = z + 1. More recently Minnhagen et al [49] 
have found a different value 

.(T)^^-l (5.24) 

by combining static and dynamic scaling arguments. They then link this expression to a 
temperature dependent dynamical exponent 

z{T) = a{T) - 1 (5.25) 

Interestingly enough, both expressions yield the same value a = 3 at Tbkt, since (see section 
1): 

J{Tbkt) = '-^^^ (5.26) 

vr 

Bormann [70] has established a link between (5.23) and a more involved recombination 
scenario - a collective <Cpartner transfer^ between bound pairs - which is supposed to take 
place when the vortex density is not too small. A recent analysis of I-V data obtained for 
various 2D superconductors [71] has produced a much larger exponent on the order of 5 or 
greater. Such a value might be justified [35] by introducing a <Cthermal escape time3> for 
a VA pair out of their logarithmic interaction potential into the generalized friction function 
determining the dynamic vortex correlator (see section 5.3). On the other hand it has been 
pointed out [51,94] that finite size effects play in important role in determining a{T). 

Shenoy [16] has developed a dynamic renormalization group approach for J J As. The eval- 
uation of the dynamic current-current correlation function Fc{lv), defined in (4.16), is based 
on the Fokker-Planck equation (4.1) for the probability distribution of the phases. Trans- 
forming the phase action (Hamiltonian (1.2) plus friction term) to vortex variables allows 
borrowing results from static BKT-scaling and imposing it onto dynamics, by introducing 
an extra frequency scaling. The dynamic conductivity cr(w) of the array is then expressed 
as an integral over length scales I of a scale dependent conductance a {to, I). This approach 
was somewhat refined later on [107]. While the aproximations for a{uj, I) in the first calcula- 
tion [16] essentially yielded a Drude-like frequency dependence, a more refined treatment of 



allowed to reproduce the MP form (5.6) in some intermediate frequency range, and thus to 
understand this anomalous behaviour from a more microscopic point of view. 



6 Disorder effects 

6.1 Linear dynamics of disordered arrays 

In the considerations of section 5 the energy barrier that a moving vortex has to cross when 
moving from one plaquette of the array to a neighboring one (see section 3.1) has been 
neglected. This approximation must be abandoned when the array deviates from perfect 
translational invariance. Different types of disorder can be considered. 

a. For bond disorder the parameters characterizing the individual junctions, namely 
the Josphson couplings Jw, the junction resistances Rw and the mutual capacities Cw, are 
random. This type of disorder is realized when the positions of the individual superconduct- 
ing islands deviate in a random way from their ideal lattice positions or if they have random 
size, since Jw, Rw and Cw depend on the distance between the adjacent superconductors 
and on their geometry. 

b. For site disorder the properties of the superconducting islands, ie their resistance 
and capacitance to ground, are random. More simply, in a ^ dilute ^ array certain islands 
are totally missing. 

Dynamic measurements have in particular been done on dilute, or <C percolative arrays 
in which only a fraction p of superconducting islands is present [95]. The experiments were 
done for p- values somewhat larger than the percolation threshold pc = 0.5 where the disorder 
has a marked influence on the properties of the array. Indeed for p < pc the average array 
has no more connected path joining opposite boundaries and superconductivity is completely 
suppressed. At temperatures that are low compared to the Tbkt (which in a dilute array is 
lower than for the regular counterpart, due to the absence of certain bonds, see, for example, 
section 17 of reference [1]) the equations of motion can be linearized, taking into account 
only small amplitude <Cspin wave^ excitations. The problem at hand is then the same as 
determining the vibrational modes of a disordered harmonic solid, which is a well studied 
field [72]. In reference [73] the impedance Z{ijo) of such a dilute J J A has been determined 
by evaluating the dynamic voltage correlation function (4.20) by different methods. The 
coherent potential approximation [CPA) has proven to be a useful approach for treating 
dynamic disorder problems. On the other hand Mori's method, applied in section 5.3 to 
vortex dynamics, is also useful. The inverse of the impedance is expressed by an effective 
coupling function K{uJ,p) depending on frequency and on the dilution fraction p : 



corresponding to the form (4.41) or to the effective inductance L and resistance R of (4.18). 





As a main result these calculations reveal the existence of three different frequency 
regimes: 



a) For very low frequency inductance and resistance are frequency independent. Such low 
frequencies effectively show the long wavelength response of the dilute array. The details of 
the disordered structure are averaged out and the <Cmacroscopic^ behaviour which is seen 
corresponds to a homogeneous superconductor. However, 1/L and R no more correspond to 
single junction parameters. They are reduced and go to zero when the percolation threshold 
Pc is reached: 

(6.2) 

R{0,p) oc (p-pcY 

where t is the critical conductivity exponent introduced in percolation theory [74]. For the 
2D lattice at hand t = 1.3, whereas the more simple <C mean field ^ treatments like CPA 
yield t = 1[73]. 

b) In an intermediate frequency range the dynamic conductance varies like a power of 
frequency : 

G{uj) - c^"" (6.3) 

reflecting dynamic scaling that results from the self-similar structure of the array near the 
percolation threshold. The precise value of u depends on the type of disorder and on the 
method (such as CPA or Mori) used for the calculation. CPA applied to bond disorder and 
bond dissipation yields n = 0.5, which corresponds to the exact result found in percolation 
theory [75] by taking into account the self-duality of the random network complex impedances 
by which the dilute J, J A can be modelled. On the other hand site dissipation yields frequency 
exponents for 1/L and R which differ from each other and are slightly smaller, respectively 
larger, than 0.5, their sum being, however, equal to 1. The cross-over frequency where regime 
a goes over into regime b is itself a function of the dilution parameter p : 

u,{p) ^{p- pcf' (6.4) 

This regime corrsponds to excitations which have a wave length smaller than the perco- 
lation correlation length and are therefore a direct witness of the disorder in the array. In 
the field of disordered harmonic solids such modes are called ^ fractons 

c) At even higher frequencies 1/L and R are again frequency independent and take the 
original values of a single junction multiplied by p. The dynamic response is now governed by 
excitations of single bonds and the relevant parameters correspond to averages over the array. 

Experimental data obtained on dilute arrays will be presented in section 8. 
6.2 Vortex dynamics of disordered arrays 

Deviations from a regular lattice structure will have an influence on the vortices of an array 
through the Peierls force or pinning potential in equation (3.5), which has to be reintroduced 
in the equation of motion (3.6) or (3.9) for the vortex centers. The equilibrium arrangement 
of the (thermal or field induced) vortices and antivortices will correspond to a minimum of 
the free energy in the given random pinning potential landscape. It is generally accepted that 
weak disorder results in some sort of glassy vortex state [76,77]. At sufficiently low tempera- 
tures such a vortex configuration would then be trapped by the underlying potential, and the 
system should have a vanishing linear resistance, ie it would be truly superconducting. This 



should be valid for dilute JJAs, as they have been described in the preceding subsection, 
as well as for ultrathin superconducting films [78] that are intrinsically disordered. There 
is, however, an important difference between 2D and 3D systems : in two dimensions the 
V-system has thermally excited defects, in particular dislocation pairs, which allow for local 
flux flow (and thus ohmic resistance), even though any motion of the V-A-structure as a 
whole is quenched by the pinning potential. Such a motion can, however, be very slow, so 
that on the time scale of a given experiment the V-system still appears to be frozen. Recent 
measurements showing such a <C dynamic freezingS> will be presented in section 8 [78]. 

In the percolative arrays described in the preceding subsection the sign of disorder has 
also been detected in the vortex dynamics. Above the BKT transition the vortex contribution 
to the dynamic impedance Z{u>) shows power law behaviour with a temperature dependent 
exponent. These data [98] will be interpreted at the end of subsection 8.2. 

7 Theoretical description of experiments probing equilibrium 
dynamics 

Three major experimental techniques have been developed in order to get insight into the 
dynamic behaviour of JJA's at or near equilibrium: 

a) I-V characteristics (V-response to static driving fields) 

b) Flux noise (spontaneous V fluctuations) 

c) Two-coil induction (V-response to weak dynamic driving fields) 

Before reporting on experimental findings we briefly explain the link between the measured 
quantities and the theoretical techniques exposed in section 5. 

7.1 Current-voltage characteristics 

Below the BKT transition there are no free vortices and the dc resistance is zero. However, 
an applied current will act on the vortices according to equation (3.6) which leads to a finite 
probability for a pair to unbind and thus to dissipation. Experimentally this phenomenon 
can be observed through the non-linear current-voltage characteristics 

V oc P^^^ (7.1) 

The form of the temperature dependent exponent a{T) for T < Tbkt can be found setting 
the escape rate over the energy barrier for a pair equal to the probability for recombination. 
Details can be found in the original work of Ambegaokar et al [34] or in the review by Newrock 
et al [1]. The result was given in (5.23). Thus, for T < Tbkt, the I-V relation is non-ohmic, 
while above the unbinding temperature the presence of free vortices leads to ohmic dissipation. 

As mentioned in section 5.4 dynamic scaling leads to result (5.24) which is different except 
at T = Tbkt- Bormann [70] has refined and extended the picture of this balance between 
pair breaking and recombination. At not too low densities of vortices charge transport is 
supposed to occur predominantly through direct particle transfer from one bound pair to a 



neighboring one. Taking into account such <^ partner transfer ^ Bormann indeed obtains 
the I-V exponent of [49] thus providing a more physical insight than simply dynamic scaling. 
However, in contrast to reference [49] Bormann's result for the exponent a(T) goes over into 
the <C classial ^ form proposed in [34] for very low currents. Thus the apparent contradic- 
tion between the two forms of [34] and [49] for a{T) is resolved. 

Theoretical values for a(T) resulting from numerical simulations will be presented in 
section 8 and experimental findings for the I-V characteristics of J J As will be shown in 
section 9. For an adequate comparison between theory and experiment one has to bear 
in mind that finite size effects can influence the BKT scenario in particular through the 
presence of some free vortices which are not the result of BKT unbinding, but rather of the 
finiteness of the system. Such finite size effects have been discussed by different authors, see 
[51], [71], [79]. An observable effect of this phenomenon will be discussed in section 8.1. 

7.2 Two-coil induction 

These measurements use a drive coil of turns of radius R£, and a pair of astatically wound 
receive coils of turns of a smaller radius Rf( mounted at some small distance d from the 
plane of the array (see figure 2). In order to study the dynamic response of the array an 
ac magnetic field is created by the drive coil, which excites the vortex system. The resulting 
motion of the vortices creates a voltage signal 6V in the receive coil. This technique was 
pioneered by Hebard and Fiory [80] for superconducting films and has proven to be very 
successful for studying the dynamic aspects of the BKT transition [1,81-85]. 



Figure 2 :Two-coil set-up used in the measurement of the dynamic conductance of the 
array, described in section 7.2. 



The measured signal 6V can be related to the dynamic conductance of the array by a 
straightforward calculation. Since the latter can be found in various references [81,83], we 
only sketch the main steps without giving an explicit derivation : 

- An ac current Id{^^) of frequency u sent through the drive coil produces a time dependent 
vector potential A{uj) 

- The time derivative of A, in the plane of the array, is an electric field E(u;) acting on the 
latter 

- An electric current density = G{uj)E{u!) is induced in the array that acts on the vortices 

- The resulting time dependent flux threads through the receive coil and can be detected by 
the induced voltage 6V{uj). 

Thus the signal 6V{ij) is linearly related to the drive current Id- 

6Viu;) = iulDiu;) r dx—^^^ (7.2) 

Here h is the sum of the closest distances of the drive and the receive coils, respectively, 
from the array plane and the function M{x) takes into account the geometry of the experi- 
mental set-up 

M{x)=7Tfio^^M—x)Ji{—x)e ^ _ (7.3) 

SHd, ShR being, respectively, the distance between the turns of the drive and the receive coils, 
in units of h and Ji is the first order Bessel function. In order to bring in vortex dynamics 
into the expression (7.2) one then has to use a link between the electric conductance G of 
the array and the vortex dielectric function e(aj), such as (4.45). More details concerning the 
derivation of (7.2) and (7.3) can be found, for example, in [83]. 

Even though finite frequencies are involved this technique allows to gain interesting exper- 
imental insight into the i^A'T-transition. Moreover it has been used to detect the frequency 
dependence of e{uj). Results will be presented in section 9. 



7.3 Flux noise 

Currents fiowing in the array produce a magnetic field that extends into the space above 
and below the plane of the array. This field can be detected by measuring the corresponding 
magnetic fiux through a given surface. In the usual experiment a circular wire of radius R 
delimiting such a closed surface is fixed at a distance d from the J J A. If the field, and thus 
the flux, varies with time a voltage is induced in the wire according to Faraday's law. In the 
absence of an external excitation (which was used in the two-coil experiment) the time varia- 
tion of the flux is due to thermal fluctuations in the array. Such a ^flux noised experiment 
measures 5*0 (w), the frequency Fourier transform of S(j,{t), which is the dynamic correlation 
function ((/)(f)(/!)(0)) of the time dependent flux <p{t) through the area surrounded by the wire. 



The usual theoretical analysis of flux noise attributes the fluctuating currents in the array 
to moving vortices. The flux (p is expressed as an integral over the dynamic charge density 
correlator of the vortex system in the plane by means of the following steps (we give some 
more details here since the relevant references do not give a complete derivation): 



is the line integral of the vector potential A along the wire 

ds ■ A (7.4) 



- The vector potential A which is produced by the electric current density j flowing in 
the array is found by solving Maxwell's equation 



A(r) = ^toj (i2A,|,e*k|,r,|£_^K,„(k|,) (7.5) 



Here the position vector r = (r|j,r_|_) is decomposed into its components in the plane of 
the array and perpendicular to the latter and is the 2-dimensional Fourier transform of 
the areal current density in the array. 

- For a single vortex the current K, circulating around its center, is given by 

K(r||) = |j(Ve(r|,)-|A(r||)) (7.6) 

Equations (7.5) and (7.6) have to be taken into account simultaneously which leads to 
magnetic screening in the final expression for the vector potential: 



" A||/cii 



The effective (in-plane) screening length Ay is given by 

A. 



2e2j/io 



- Assuming that the relevant contributions to the phase gradient are due to moving vortices 
the current density K is expressed by the Fourier transform of the density of the latter : 

Ki„(k||) = -27ri(^)J^^^p,(k||) (7.8) 

- Introducing expressions (7.5) to (7.8) into (7.4) leads to an explicit expression of the 
flux (j) through the plane of the wire in terms of the vortex density in the plane of the array : 

= ^i2^o(-27ri)^J r dkMik.R)e-h'i^^^P^ (7.9) 
n Jo 1 + A||A;|| 

where Ji is the first order Bessel function bringing in the geometry of the pick up wire. 



- The flux noise correlation is given by 

S^u;) = Jdt e'^\mm) (7.10) 

- The statistical average concerns the vortex density showing up in statement (7.9). It 
leads to the dynamic density correlator of the V-A system introduced in equation (4.24). The 
latter is diagonal in the wave vector k (as it is defined in (4.24)). Thus the S^(lo) is expressed 



by the Fourier transform of the dynamic vortex correlation function which, in turn, is related 
to the Laplace transform (4.24) by 



J dt e'^\p,ik,t)plik,0)) = 2Re^pp{k,-iu;) (7.11) 
The final result for the flux noise spectrum reads 



S^{lo) = {2tt'Rp,Y{-Y{J\Y dk l^^^^^^j^y 2Re^pp{k,-iu:) (7.12) 

In case several concentric wires are set up on top of each other (7.12) has to be replaced 
by a sum of terms of the same structure, but each with the corresponding distance d„ from 
the array, showing up in the exponent. Results obtained through expression (7.10) by using 
different forms for the dynamic vortex correlation function (j)pp{p, z) will be shown in section 
8. They will be compared with experimental data in section 9. 

In order to build a bridge between flux noise and electrical conductance one can use the 
links between charge density correlator, dynamic susceptibility and dielectric function shown 
in (4.26,27) in order to rewrite (7.12) as 

2d ^2f'^e2fj..2f'^r,Ji{f^R?e-^'"^2kBTk'^ 1 



S^c) = {27.'RpoY{-Y{JhY / dk I) ' --f^ — (7.13) 

n " Jo k(1 + A||A;)^ uj 2T:q^n e^iu}) 

When the pick-up coil is very close to the plane of the array {d small) only small /c-values 
contribute to (7.11) which can be approximated by 

S^{uj) = C-Im{-^) = C'ReGsiio) (7.14) 

relating thus - through (4.45) - flux noise approximately to the long wave length conductivity, 
or - more precisely - to its superconducting part [48] . 

The above flux noise calculations have been based on taking into account only the super- 
current in the array (see equation (7.6)) and by relating the phase gradient in (7.6) directly 
to vortex configurations. The link between flux noise and conductance can be made more 
general [86] by expressing first the flux correlator (7.10) by the correlator of the total array 
current and by relating the latter to the array conductance via the considerations exposed in 
section 4 on linear response. 



8 Experimental data 

We present experimental data for the three types of <Cstandard^ experiments for which we 
have related the measured observables to theoretically accessible quantities in the preceding 
section. 



8.1 Current- voltage characteristics 

Various measurements of the relationship between an applied current and the resulting volt- 
age drop, yielding the exponent a(T) in (7.1), are cited in the review by Newrock et al [1]. 
There (in section 11) the interplay between different relevant length scales implied in such 



an experiment, where currents of different intensity are used, is also explained in detail. The 
jump from the non-ohmic behaviour given by a value of a(T) that differs from 1 to the high 
temperature linear I-V relation is clearly bourne out by experiment. 

As explained in sections 4 and 6 a new T-dependence [49] of the coefficient a, differ- 
ing from the <C classical ^ form has recently been proposed. Weber et al [64] have re- 
analysed I-V data (see Figure 3) obtained on Hg-Xe alloy films [88], Bi2Sr2CaCu20x[89] 
and BiisPboASr2Ca2Cu30x [90] single-crystal films, finding indeed good agreement with the 
new form of a(T) of Ref. [49] and [70]. Herbert et al [79] have come to the same conclusion 
for proximity-coupled and high- Tc weak-link arrays. They also find evidence for the finite 
size effects mentioned at the end of section 7.1. Indeed, the I-V curves for very weak currents 
bend over into a linear form even below the transition temperature, since the free vortices 
present in a finite array, lead to ohmic behaviour. 



Figure 3 : Exponent a{T) of the I-V characteristics (see equation (7.1)), according to 
Ref. 87, plotted as a function of X = T/Tc- The different curves have the following meaning 

- Plusses : Reference 88 

- Other symbols : Monte Carlo results for an array of length of 16 (stars), 24 (open circles), 
32 (filled circles) and 48 (triangles) unit cells. 

The dashed line corresponds to the universal value a = 2, at the critical temperature. 
The inset shows experimental data for Hg—Xe alloy films, Ref. 88 (plusses), Bi2Sr2CaCu20x 
single-crystal films, Ref. 89 (filled squares) and Bii,QPhQ,iSr2Ca2Cu^Ox single-crystal films, 
ref. 90 (open squares). 



8.2 Two-coil induction 



The two-coil mutual induction technique, from which the dynamic conductance G{lv) of the 
array can be extracted (see section 7), has extensively been used in order to study the BKT 



transition. One of the basic features is the temperature dependence of the phase stiffness or 
the hehcity modulus A, reaching a universal value at Tbkt before jumping to zero. As ex- 
plained in section 4 the low frequency limit of G{ijj) is directly related to A. Ref. [81] presents 
data obtained on a proximity-effect J J A consisting of Pb-islands on a Cu substrate. The 
voltage induced in the pick-up coil is analysed according to section 7.2 and the correspond- 
ing dynamic dielectric function e(a;) is extracted for different temperatures in the vicinity 
of Tbkt- The fall-off of Re{e~^) near Tbkt is characteristic of the ^insulator-conductor 
transition^ of the vortex system and Im{e~^) shows that this transition is accompanied by 
dissipation. Going to higher frequency shifts the transition region to higher T, since larger to 
tests shorter length scales for which the unbinding of V-A pairs happens at higher tempera- 
ture T^j, which can be interpreted as an <C effective frequency dependent critical temperature 
In order to extract the helicity modulus A which is a thermodynamic response coeffi- 
cient one has to extrapolate to zero frequency. This can be done by identifying iie(e(a;,T)) 
for a given frequency lo with the static dielectric function e{l^,T) for a length scale cor- 
responding to oj. The T-dependence of A obtained in this way is in good agreement with 
a renormalization group calculation integrating the scaling equations up to a length of the 
order of S,{T^), the BKT correlation length at T = T^^ where pairs of length ^(T^^) begin to 
dissociate. The extrapolation to infinite length scale agrees very well with the universal jump 
of A at Tbkt prredicted by theory. More details concerning this scaling procedure used for 
getting at A can be found in [81] and [91]. 

It is important to mention that whenever experiment and theory is compared in this way 
one has to take into account the T-dependence of the Josephson coupling given, according to 
(1.3), by the one of the critical current 

h{T) = 1,(1 - ^fexp{-^) (8.1) 

see [92], where Tbcs is the BCS transition temperature of the individual islands of the array, 
separated by the distance L and ^„ is the effective coherence length of the normal metal 
forming the junction between two neighboring islands. One can then compare experimental 
data with XY-model calculations by using a dimensionless temperature parameter 

' /cfiT 2ekBT 

Moreover, spin wave excitations produce an additional renormalization of the Joseph- 
son coupling. The latter can be taken into account by using an <^ effective Coulomb gas 
temperature » [93] 

f 

T = ^ (8.3) 

The two-coil technique has also allowed to shed light on the dynamics of 2D Coulomb 
systems. In section 5 two phenomenological approaches to 6(0;) have been presented : Drude's 
model and Minnhagen's phenomenology (MP). One might expect that above Tbkt, where 
free vortices are supposed to dominate the response to a driving force, Drude's model would 
be appropriate. It thus came as a surprise when measurements on weakly frustrated triangu- 
lar J J As [82] revealed a frequency dependence of €{uj), see figure 4, which is well described 
by MP and support the idea of a non-analytic dynamic mobility, as explained in section 
5.2. Triangular arrays are better suited for studying the intrinsic dynamics of interacting 



vortices since the barriers due to lattice pinning (section 3.1) are very low [23]. One has to 
bear in mind that the system of field induced vortices in a frustrated array should actually 
be described as a <C one-component plasma in a neutralizing background ^ (OCP) rather 
than the <Cneutral two-component plasma^ (TCP) of the thermally excited V-A system in 
an unfrustrated array. Thus, theoretical approaches developed for the TCP, as in section 5, 
are not immediately applicable. However, numerical simulations of frustrated arrays (section 
9) have shown that the two types of plasmas behave the same way. 



Figure 4 : Real and imaginary parts of the inverse of the vortex dielectric function mea- 
sured in reference 82 on a JJA in a small magnetic field. The two functions are plotted as 
a function of the reciprocal frequency. The full lines are fits to the data using Minnhagen's 
expressions (5.6), whereas the Drude form is given by the dashed lines. 

It is clear that such a vortex mobility of the form (5.6) cannot hold down to arbitrarily 
low frequency above the vortex unbinding transition. Indeed, the electrical conductance G{uj) 
is related to the vortex mobility according to (4.46). The vanishing vortex mobility (5.6) for 
u; — > of MP would therefore imply that superconductivity is maintained even above Tbkt- 

Another interesting application of the two-coil technique has shed light on the dynamics 
of percolative arrays [95] discussed theoretically in section 6, see figure 5. The data confirm 
that the transition termperature goes to a lower value when the concentration of missing 
sites increases, ie when the percolation threshold is approached. Moreover, both the inverse 
inductance and the resistance show very clearly two of the three frequency regimes discussed 
in section 6 : for low frequencies both quantities have constant values given by (6.2), whereas 
above some cross-over frequency the fractal behaviour is observed. The critical exponent 
t ~ 1.4 is extracted from the data which is in good agreement with the prediction t ~ 1.3 
for percolation in two dimensions [96]. The frequency dependence in the fractal regime is 
compatible with u ~ 0.55, see equation (6.3) which is somewhat larger than what is expected 
from duality arguments for bond percolation [75] and for a mean field approximation. On the 
other hand, different types of approximations used in the theoretical description of disordered 



arrays in Ref. [73] also lead to u values that differ from 0.5. 



Figure 5 : Log-log plot of the frequency dependence of the inverse kinetic inductance L~ 
and resistance R (according to equation (4.18)) of a dilute array close to the percolative limit 
{p = 0.51) from reference 95. The dashed lines indicate uj^^"^ power laws. 

Dynamic response measurements on the dilute arrays near their perconlation threshold 
have also been extended to temperatures near and above the expected BKT temperature [98] , 
where vortex motion is expected to dominate over the linear spin wave dynamics discussed 
before. The contribution Zi,{u;) of the thermally excited V-A-system to the impedance Z{uj) 
is extracted from the measured values by representing the impedance of the array, in the spirit 
of the RSJ model, by two parallel channels - an inductive one for <^ superfluid ^ conduction 
and a dissipative Ohmic path - and by assuming that the vortex response renormalizes the 
superfluid part. A similar decomposition as (4.18) fo the dynamic conductance is used for Zy 

Z,a{uj) = Ry{uj) + iujLy{uj) (8.4) 

(whereby one has to bear in mind that resistances and inductances introduced, either by 
(4.18) or by (8.4) are not identical ). Above Tbkt both, R^ and L^,, show a power law 
behaviour with frequency 

Ry{uo) oc , L^w) oc (8.5) 

with an temperature exponent. The latter is close to 1 at Tbkt and decreases slowly when 
T increases. Such a T— dependent exponent can be obtained by assuming that the relevant 
contribution to the dynamic response comes from pairs and by calculating their polarisability 
as indicated in subsection 5.1. However, a simple calculation using statement (5.3) yields an 
exponent u{T) that decreases much faster with rising T than the ones observed [98]. Thus, 
one should probably revise the ideas of subsection 5.1 for disordered arrays. The response at 
temperature much higher than Tbkt again points to the motion of free vortices, the latter 
being activated across the pinning barriers. 



Glassy behaviour of the vortex system, mentioned in section 6.2, has been seen in mea- 
surements of the dynamic impedance Z{uj) in ultrathin Y Ba2Cu2,OT{Y BCO) films in a per- 
pendicular magnetic field [78]. At high enough temperatures the dissipative component of 
Z has an activated T-dependence resulting from activated motion of single (field induced) 
vortices out of local pinning wells. With decreasing temperature Z crosses over to a more 
or less T independent value. The cross-over temperature T*, that depends weakly on uj, 
turns out to be much higher than a reasonable estimate for the melting temperature of a 
2D V-lattice, due to unbinding of dislocation pairs. This points to the effect of pinning on 
the V-system. The frequency dependent response of such a pinned elastic manifold has been 
studied theoretically by various authors [99]. The corresponding predictions agree well with 
the measurements on YBCO films. It would, of course, be interesting to see in more refined 
experiments whether T*{uj) decreases without limit when the frequency oj becomes lower and 
lower. This would be a true sign of the fact that the 2D <C glassyness ^ is really a dynamic 
phenomenon. 

8.3 Flux noise 

The flux noise spectrum S(f,{u)) has been measured in overdamped J J As by Shaw et al [68] 
and more recently by Candia et al [100], as well as in hi^-TciY ECO) films by Festin et al 
[101]. Shaw's data cover the range from T = 1.825i^ to T = 2.379il' for a Nh on Cu array 
with Tbkt ~ 1.63-fC, see figure 6. At each T, S^{uj) is white for u < ujc and varies like 
1/u! for uj > uJc ■ The cross-over frequency uJc is proportional to the inverse square of the 
Kosterlitz-Thouless correlation length with the T dependence as given by (1.12). This result 
is consistent with dynamic scaling (see section 5.5), for a dynamic exponent z = 2 and with 
6 ~ 2 in expression (1.12). The white noise level increases for decreasing T as Cbkt{T)'^ ■ 
The authors conclude that, even at the lowest temperature shown (ie., for T as close to Tbkt 
as possible), S,bkt is on the order of the size of the pick-up coil and, as a consequence, the 
prefactor in (1.12) is actually smaller than the lattice constant of the array. The fact that 
scaling is obeyed down to such a microscopic length scale is somewhat surprising. 

The data of ref. [101] confirm previous observations [102] showing that high- Tc films also 
show typical 2D vortex fluctuations. Data are shown for 500 A thick YBCO films from 
T = 84.84i^ to T = 84.89i^, just slightly above Tc. At low frequency 5'^(i^) is again white, 
near some frequency u!o the data have a common tangent corresponding to 

S^iuj) ~ 0.-1-05 (8.6) 

before crossing over to 1/cj^/^. However, contraray to what has been seen by Shaw et al 
[68] the noise spectrum does not have a true l/w dependence over any appreciable frequency 
range. Thus, the two materials {J J A, respectively high Tc film) have a different characteristic 
behaviour : while the J, J A shows an impressive 1/u; noise, the film data only have a l/w 
tangent whereas the high frequency behavior is 1/cj^/^ . The same is found in [102] for 
Bi2Sr2CaCu20s+x films. 
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Figure 6 : Flux noise spectrum versus frequency for different temperatures of a JJA 
consisting of Nb islands on Cu (ref. 68). The dashed hues have slope -1 and 0, respectively, 
signalling white noise and noise in different frequency regimes. The inset shows the 
derivative of the current voltage characteristics, dV/dl, as a function of temperature. The 
BKT temperature of the array is at l.GSiC. 

The most recent measurements on a proximity effect Ph/Cu/ Ph triangular JJA [100] 
confirm the l/w-dependence of above Tbkt (see figure 7), but the analysis of the low 
frequency white noise level has been interpreted in terms of an exponentially activated (total) 
number n{T) of vortex excitations. This would show that the vortex mobility determining 
(7.13) through (4.28), (see also (9.5)), which is supposed to go to zero when the critical 
temperature is approached, does not really influence the white noise level. This may be a 
finite size effect (see section 9.3 where results of theoretical flux noise calculations will be 
presented). Candia et al [100] have also extended their measurements to temperatures close 
to, but below Tbkt- In the frequency range covered by the measurements the noise keeps to 
behave like but, for a given frequency, the level decreases when T goes down. 

The ubiquity of the l/w-dependence is somewhat surprising. As shown in section 9 it 
can be reproduced when pair dynamics is taken into account, but with an exponent in the 
probability distribution for the relevant pair distances in statement (5.2) which is frozen at 
T = Tbkt- Thus it seems that the array <^ stays critical even when the temperature 
differes from its critical value, as it has been remarked by Shaw et al [68]. In the context 
of disordered arrays [98], the corresponding exponent has also been extracted from dynamic 
data and shown to be much less temperature dependent than statement (5.2) would imply. 
As a matter of fact, the 1/uj dependence could quite generally be a disorder effect. Indeed, 
even small variations of the distance d between neighboring superconducting islands implies 
a relatively strong <C disorder ^> in the Josephson coupling constants, since the latter depend 
exponentially on d. In such a disordered array the free vortices that can be present at any 
temperature will have to move on a ^ random substrate The corresponding dynamic 
response would typically behave like \/lo. 
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Figure 7 : Flux noise spectrum from reference 100, measured on a Pb/Cu array as a 
function of frequency / at temperatures slightly below (a) and slightly above (b) the BKT 
transition. 



9 Theoretical results 

9.1 Dynamic critical exponents 

The dynamic exponent z, introduced and discussed in section 5.5, has been determined 
through various numerical simulations. Tiesinga et al [46] solved the equations of motion 
for the phase field using alternatively TDGL and RSJ dynamics, as explained in section 
5.4. The corresponding value of z is obtained by studying the scaling form (5.21) of their 
correlation functions, ie by identifying the link between the relevant relaxation frequency and 
the BKT correlation length. Interestingly, they find that z is different for the two types of 
dynamics : z{TDGL) ~ 2 corresponds to what is expected, whereas z{RSJ) ~ 0.9 points 
to a different type of critical dynamics(see the insets in figure 10 that will be discussed in 
connection with flux noise). This would be an example of a system with a given static criti- 
cal behaivour [BKT in this case) which goes over into different dynamic universality classes, 
depending on the precise equations of motion used to describe time evolution. 

The scaling analysis of the TDGL and RSJ simulations of Kim et al [53], who used 
periodic and fluctuating twist boundary conditions [FTBC, see section 5.4), yields z = 2 



at the transition temperature for both types of dynamics for FTBC. This value is found 
by finite size scaUng (extrapolation of calculated parameter values to the thermodynamic 
limit) of the array resistance. However, they observe that finite size scaling with PBC gave 
different results, or - more precisely - that PBC are actually inadequate for this procedure, 
since the relaxation time which reveals the exponent z is proportional to the macroscopic 
resistance that is actually zero for periodic boundary conditions. This shows once again that 
correct boundary conditions are crucial whenever transport coefficients are considered, which 
describe the response of the system to an external influence that couples to the boundary. 
Below Tbkt the value of z is slightly larger than 2, which is in agreement with the finding 
for a 2-d lattice Coulomb gas with Monte Carlo dynamics [64]. 



9.2 Vortex dielectric function 

The fact that the vortex dielectric function, measured by the two-coil technique described 
in the previous section, has a frequency dependence resembling more the <CMinnhagen 
phenomenology^ {MP) than the classical Drude behaviour, has initiated both analytical 
and numerical work aiming at understanding this <C anomalous^ dynamic response of J J As. 
Such an effect is likely to be related to vortex-vortex interaction, to pinning by the under- 
lying lattice or to the fact that vortices are not just point particles, but singular excitations 
of the phase field which change shape when they move and when they interact with each other. 

Coupling to spin waves [103,104], which are the <CBosons^ mediating the Coulomb inter- 
action between vortices in the framework of the effective electrodynamics described in section 
4.2, has been shown to lead to a vortex dielectric function of the form given in section 5.2. 
The non-analytic low frequency dependence, which is characteristic of MP, is an <C infrared 
divergence ^ that can be understood by counting powers of wave number in the relevant 
integrals. It is the combined effect of the spin wave propagator and the dimensionality of the 
array. Other authors [105,106] have also found non-analytic effects in vortex dynamics, for 
instance in the vortex mass of arrays with capacitive effects, caused by coupling to spin waves. 

Capezzali et al [107] have treated the problem of interacting vortices, both by a refined 
dynamic scaling (a generalization of the earlier approach in ref. [16] which did not yield any 
MP features) and by a Mori calculation of the vortex density correlator (see section 5). They 
find three different regimes corresponding to different ranges of the scaling variable 

y = A- (9-1) 

z being the dynmic exponent discussed in section 5.5 and ^ the correlation length. The ex- 
perimental data of ref. [82] fall in the intermediate regime where MP is valid. The high 
frequency regime corresponds to <C critical dynamics ^ as it had already been exhibited in 
[16]. Contrary to the ^ pure MP ^ form there is a low frequency regime where normal 
Drude behaviour should be restituted. This is an important complement to MP, since taking 
the MP expressions (5.5) to (5.7) straight down to a; = 0, and inserting them into expression 
(4.45) for the electrical conductance of the array would yield superconducting low frequency 
behaviour even above the transition temperature. 



The calculation of the dynamic charge correlator (j)pp{see section 4.2) by Mori's technique 
has more recently been generalized [108] by replacing the Coulomb potential by its screened 



form, as it is usually done in calculations for quantum Coulomb systems [40,41]. Since screen- 
ing is due to free vortices the relevant screening length (5.13) is the BKT correlation length 
(1.12) which is directly related to the inverse density of free vortices. This allows to im- 
pose static critical behaviour on dynamics. The resulting friction function, determining (ppp 
and the dielectric function has indeed a logarithmic frequency dependence, leading to MP 
behaviour, over a frequency interval that is increasing in length when the transition temper- 
ature is reached at which ^bkt diverges. This result shows that dynamics is <C normal ^ 
as long as the long range Coulomb interaction is well screened by the unpaired particles and 
it gets more and more MP- like when screening becomes less efficient (see figure 8). This 
manifests itself in a friction function that increases more and more for low frequency, making 
the motion more <C sluggish At the transition temperature this <^ critical slowing down 
3> is completed : the zero frequency friction function diverges. 
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Figure 8 : Real part of the inverse of the vortex dielectric function, as a function of 
frequency, calculated in the framework of Mori's approach taking into account the motion of 
unbound vortices moving in a screened Coulomb potential, for different temperatures above 
Tbkt{T /Tbkt = 1-02, 1.09 and 1.37 from top to bottom respectively). The frequencies are 
given in units of a scale frequency u>a = ksT/To?, involving the lattice constant a and the 
bare damping T (see equation (5.10)). The slopes of the two straight lines correspond to 
l/e(cj) ~ u) (Minnhagen phenomenology) and l/e(w) ~ u;^ (Drude behaviour). 

While these analytic calculations treat the motion of ^ free ^ vortices in the force field 
of the Coulomb gas, the calculation of the vortex dielectric function based on the pair re- 
laxation picture [34] has also been refined [35] in order to study possible deviations from a 
simple Drude behaviour. Using a length scale dependent static dielectric constant, obtained 
by numerical integration of Minnhagen's version [17] of the BKT scaling relations below and 
above Tbkt-, yields an e(u;) that is indeed closer to MP than to Drude in some intermediate 
frequency range. 



Minnhagen's group has studied JJA dynamics in great detail [55,57,58,60,63] using dy- 
namic simulations some aspects of which have been discussed in section 5.4. These authors 
have used the different description (a) to (d) of J J As quoted at the beginning of section 5.4, 
and they have treated both unfrustrated and frustrated arrays, in which field induced vortices 
are also present. They also replaced the simple cosine interaction of the phase Hamiltonian 
(1.2) by the form 

H = J:J[1-cos''\^)] (9.2) 

This coupling reduces to (1.2) for p = 1, but for larger value of p it approaches a <C 
potential well ^ favoring the creation of vortices [109]. The dielectric function obtained by 
calculating the relevant phase correlation function indeed follows closely the MP prediction 
(see figure 9): contrary to the simple Drude form Re{l/e) varies linearly with oj, Im{l/e) 
has a long <C flat ^ region, and the peak ratio (5.8) is close to tt/2. Deviations from MP for 
small uj, leading back to Drude, appear for p > 1. Frustrated systems [54,57] show similar 
behaviour : when the magnetic field (and with it the effective array resistance) is increased 
at a fixed temperature the response crosses over from normal to anomalous. In ref. [57] it 
was observed that increasing the magnetic field - ie the frustration of the array -not only 
leads to an obvious higher density of field induced vortices, but it also increases the number 
of V-A pairs. This seems to corroborate the idea that V-A correlations are responsible for 
the anomlous dynamics observed in simulations and in experiments. The screened Coulomb 
potential used in the Mori calculations [108] may indeed be one of the manifestations of the 
presence of a sufficient number of thermal vortex excitations. 



Figure 9 : Real and imaginary parts of the inverse vortex dielectric function of an un- 
frustrated array using an interaction potential with p = 2 in equation (9.2), The open circles 
(iie(l/e)) and the filled circles (/m(l/e)) are the results of dynamic simulations using TDGL 



dynamics (method b of section 5.4) of Ref. 55. The dashed and full curves are fits to the 
numerical data using the MP relations (5.6). 

The effect of (regular or random) pinning on the motion of field induced vortices of 
weakly frustrated JJAs has been investigated by Monte Carlo simulation [65] of the 2D 
one-component lattice Coulomb gas (OCCG). Anomalous MP behaviour of the dynamic 
dielectric function (see subsection 5.2) is found, and it seems to be related to pinning. It is 
well developed when sufficiently strong random pinning is present. In the absence of random 
pinning it is restricted to temperatures close to the depinning of the vortex system from the 
underlying lattice structure. On the other hand, the continuum OCCG, which experiences 
no pinning at all, shows no anomalous response. In contrast to this, the neutral 2D two- 
component Coulomb gas [62], which respresents the VA-system in an unfrustrated array 
shows MP behaviour below the BKT transition. One may therefore wonder whether in the 
latter case vortices of one sign of charge act as pinning centers for antivortices, and vice- versa 
[62]. 

9.3 Flux noise 

Dynamic simulations have also been used for evaluating the flux noise spectrum either by 
evaluating the full integral (7.12) or by using approximate forms of the latter, such as (7.14). 
At low frequency the expected white noise is found with a level that goes up when T ap- 
proaches Tbkt from above. The behaviour for higher frequencies varies between the different 
authors. Some calculations [56], [110], simulating TDGL dynamics for a 2D XY model with 
a phase coupling of the type (9.2) yield 

S4,{oj) ~ (9 3) 

for some intermediate frequency range and S(j,{u)) ~ w"^ for high frequencies. The broken 
power of (9.3) can be understood as the sign of conventional random walk of vortices across 
the boundary of the pick-up coil [110,111]. When the curves bend over from white noise to 
the power law (9.3) a common tangent can be drawn to the curves which has a slope close 
to However, there is no extended regime showing the measured noise (see section 

8.3) . 

Other dynamic simulations [46] find an intermediate frequency range, above the white 
noise regime, where the spectrum follows in power law : (l/cj)" (see figure 10). The exponent 
a is T-dependent and differs slightly between different types of dynamics used (see section 

5.4) : 

0.85 < a < 0.95 for TDGL dynamics 

(9.4) 

1.17 < a < 1.27 for RSJ dynamics 

which, again points to different 'critical dynamics' where different equations of motion are 
used. Both cases are relatively close to the measured dependence, but - as the authors 
state - the TDGL dynamics is closer to the measured l/co dependence. This is somewhat 
surprising, since the <Cusual^ picture one has of a Josephson juncion and of an array is 
precisely the RSJ effective circuit model. It is possible that the arrays of Shaw et al [68], 
consisting of large Nb islands deposited on a Cu substrate have very small resistances to 
ground, such that the local dissipation term in the TDGL equation dominates the bond term 
in RSJ. Other i^S'J-simulations [112] also yield the behaviour of (9.3), both for / = and 



for frustrated arrays. The same holds for the solution of the Fokker-Planck equation for pairs 
[43], discussed at the end of section 5.3. On the other hand a calculation [113] of the flux 
noise spectrum starting from a Euclidean action for the XY-model with local damping also 
finds a cross-over from white to l/w noise. 



Figure 10 : Flux noise spectra versus frequency obtained in Ref. 112 for TDGL (a) and 
RSJ (c) dynamics for different temperatures above Tbkt- The insets of the panels (a) and 
(c) show the exponent a (see our expresssion (9.4)) as a function of temperature. Panels (b) 
and (d) show a scaling plot of the noise spectrum as a function of time. The insets present 
the dependence of the characteristic time t on the correlation length x, yielding the values 
z = 2, respectively 0.9 of the dynamic exponent introduced in subsection 9.1. 

The results obtained by the Mori approach to the dielectric function (section 9.2) have 
also been used for evaluating flux noise spectra. Preliminary calculations [108] show that the 
friction function, obtained by taking into account motion of unbound particles in a screened 
Coulomb potential, although reproducing nicely the MP behaviour in e{uj), does not lead to 
l/u> noise. These calculations are thus in good agreement with the above mentioned numer- 
ical simulations [56], [110]. Another approach starts from the averaged pair susceptibility 
presented in section 5.1, integrating the single pair susceptibility x over a pair distribution 
of the form (5.2) from the lattice constant (the smallest possible pair distance) up to the 
BKT correlation length ^ and going from x to the density correlator (f)pp with the help of 
(4.303) yields a flux noise spectrum with a white noise level and a subsequent power law with 
a temperature dependent exponent [97], similarly as in [46]. In fact, the l/w dependence is 
found for T = Tbkt- The white noise level and the cross-over frequency between the latter 
and the power-law regime varies like Taking the limit of zero frequency in the integral 
(7.10) indeed yields the following approximate statement for the white noise level 



(9.5) 



where the constant is the result of the wave number integral involved in (7.10), n.u{T) is the 
total number of vortices and 7(0) the zero frequency limit of the friction function. The data 
of reference [100] for the white noise level can be understood by limiting the integral over the 
pair distribution by the sample size L, whenever ^ > L. This explains why no more critical 
behaviour is seen in the white noise data when the critical temperature is approached : the 
relevant length is fixed by L and the only major T-dependence left in (9.5) comes from n^{T), 
which is indeed thermally activated as seen in the data of ref. [100]. 

10 Driven arrays : dynamics far from equilibrium 

J J As are interesting examples of systems with non-linear equations of motion (EM). The 
relevant EM for the phases, depending on the type of dynamics, have been found in section 
2.1, whereas section 3 ends up with an equation of motion for the vortex positions. While, in 
sections 5 to 7, we have mainly been interested in equilibrium dynamics and in response to 
weak driving fields we now turn to work dealing with the dynamic behaviour of J J As under 
moderate and strong driving currents. 

The interplay between vortex pinning due to the discrete underlying lattice structure (the 
<C egg cartoon potential ») and V-V interaction can be tested by changing the applied per- 
pendicular magnetic field, i e the frustration. In the ground state the field induced vortices 
tend to form a regular arrangement minimizing their interaction as well as the pinning en- 
ergy. The critical current needed to drive this vortex lattice across the sample is the higher 
the better the <C matching between array and vortex structure, which is the case for 
f = p/q with small integers p and q. Marconi and Dominguez have studied the depinning 
and melting of driven moving vortex lattices for small frustration [114], / = 1/25, and for 
<^ full ^ frustration [115], / = 1/2, by solving the RSJ Langevin equations of motion for 
the phases in presence of the applied vector potential and a random force which maintains a 
fixed finite temperature. The results can be summarized in a non-equlibrium phase diagram 
in the plane of temperature T^ep and drive current /. For very low current it reveals the 
equilibrium depinning transition temperature above which the vortex lattice floates freely on 
the arrays, thanks to its thermal energy. Stronger currents shift this depinning transition to 
lower T and produce phases of different properties above T^ep- At low T the moving lattice 
is still <Ctransversely pinned^ : it has anisotropic Bragg peaks, quasi-long range order and 
maintains superconducting coherence (and thus a finite critical current) in the direction per- 
pendicular to the direction of the strong drive current (one has to bear in mind that a vortex 
moves in the direction perpendicular to the drive current, see section 3.1). At higher T the 
vortex lattice is simply <C floating » . 

The presence of random pinning adds new complexity. As mentioned in section 6.2, 
an irregular underlying lattice structure forces the vortices to <Cadapt^ their equilibrium 
structure in order to minimize their free energy. In order to drive the vortex arrangement 
across the <Cpinning landscaped the applied current /, which acts as a transverse force on 
the vortices, has to exceed some threshold value Ic- Once the vortices are depinned, which 
is seen by the appearance of a finite voltage and thus a finite resistance of the array in the 
current direction, their motion will, of course, be influenced by the random pinning potential. 
Dominguez [116] has observed, in his numerical simulations, a dynamic transition occuring 
for a current Ip{> Ic)- For currents with a strength between Ic and Ip the moving ^vortex 
liquids has an isotropic structure and the turbulent flow regime is plastic. Above Ip there 



is a regime with homogeneous ^laminar^ flow in the direction of motion and with a struc- 
ture with anisotropic short range order of smectic type characterized by the appearance of 
corresponding Bragg peaks. 

Various other work has been pubhshed about the possibhty of such a dynamic regime in 
superconducting films (for references, see for example [117]). Two dynamical phase transi- 
tions for thin films are found theoretically in [118] : for a given strength of the exciting force 
the system goes again over from plastic flow, where the vortices move in an intricate network 
of channels, to smectic flow. This transition is characterized by a maximum in the differential 
resistance (the derivative of the averaged vortex velocity with respect to the drive current 
strength). At a higher driving force there is a sharp transition to a frozen transverse solid. 
Here the motion of the vortices transverse to the flow direction, which consisted of jumps 
between different channels in the smectic regime, is strongly reduced. More details about the 
temporal behaviour of the system (noise spectrum, diffusion exponents etc.) can be found in 
[118]. Moving anisotropic vortex structures have also been observed experimentally [117] in 
small high quality single crystals of NbSe^. 

Non- linear dynamic systems, such as our J, J As, can also show various phenomena related 
to the occurrence of chaos. Indeed, already a single, capacitively and resistively shunted junc- 
tion with an applied voltage or current corresponds to a damped driven pendulum, which is 
known to show very rich structures in its dynamic phase diagram. Regimes of periodic and 
quasiperiodic motion alternate with chaotic regions when parameters, such as the drive force, 
are varied. Josephson junctions coupled in series have a spatial variable in addition to time. 
They thus allow for spatio-temporal chaos. Various such systems have been studied [9] nu- 
merically. They show coherent, ordered, partially ordered, fully turbulent and quasiperiodic 
phases. When time is discretized in the coupled non-linear equations of motion for the JJA 
the latter are of the same structure as the equations of <Cglobally coupled maps» [9] . Their 
mathematical properties (transition from regular to chaotic motion, structure of attractors, 
etc.) have widely been studied. Chaotic motion in J J As should in principle be observable 
experimentally through the characteristic features of the current-voltage characteristics of 
strongly driven arrays, but not much systematic analysis seems to have been done up to now. 



11 Miscellaneous 

Josephson junction arrays are an interesting example of articifial systems. They allow us to 
profit from the microscopic phenomenon of superconductivity in order to observe challenging 
collective effects. This article is mainly devoted to the dynamic manifestation of the latter. 
The interested reader should find elsewhere more information about the static and thermo- 
dynamic behaviour of J J As in order to fully appreciate the rich spectrum of time dependent 
phenomena discussed here. For obvious reasons we had to make a choice concerning the 
topics to be presented (and this choice is certainly biased by our own experience in the field 
!). In fact, quite a variety of other aspects of JJA dynamics have not been covered. Let us 
just mention a few of them: 

Hall effect 



In addition to the ^Lorentz forceS>, mentioned in section 3, a vortex is also subject to 



a <SMagnus forced [25] when an external charge introduces a potential between the array 
and the ground. This force is perpendicular to the vortex velocity. The motion resulting 
from the combination of Lorentz force (showing up in equation (3.6)) and Magnus force leads 
to a <C vortex Hall effects that has been observed for classical and quantum arrays [25]. 
Since the Hall angle is inversely proportional to the damping constant F showing up in the 
vortex equation of motion, Hall measurements give information about the type of dissipation 
prevailing in the array. As a function of frustration, ie of the applied magnetic field, the Hall 
resistance shows a rich structure [120]. Reversal of the sign of the Hall constant, which is a 
sign that the relevant type of charge carriers changes, appears at several frustrations. 

Ballistic vortex motion 

When the mass term in the equation of motion (3.6) dominates over the friction vor- 
tices can in principle move almost freely across the sample. Such <Cballistic^> vortex motion 
has been observed, both in one- and two-dimensional arrays [1,25], in a frequency window, 
bounded from below by the presence of pinning forces and from above by various (non-linear) 
damping mechanisms setting in at high velocities. It is interesting to note that, from a more 
elaborate theoretical point of view, <Cideal^ ballistic motion is not really likely to be ob- 
servable [1]. Simulations on finite systems have to cope with boundary effects (a vortex being 
attracted by its image when it approaches the boundary). Moreover, non- linear viscosity 
can lead to localized oscillations rather than propagation of the vortex. Thus ballistic vortex 
motion definitely needs to be studied in more detail [1]. 

Shapiro steps 

In section 4.1 we have introduced a weak applied external electric field into the equations 
of motion and we have calculated the linear response of the array. For an applied voltage 
of finite strength one can use the Josephson relation (4.19) between a voltage drop across a 
junction and the corresponding time dependence of the phase difference. When the applied 
voltage is of the form 

V{t) = Vo + Vi cos uJit (11.1) 

the corresponding phase difference is easily found by integrating (4.19). Substituting it into 
(1.6) yields an statement for the resulting Josephson current for a single junction. Expressing 
the sin of a trigonometric function by a series of Bessel functions [1] shows that one should 
observe a dc super current whenever 

Vo = Vn = (11.2) 

These voltage steps are called <CShapiro steps^ [121] of a single junction. Including the 
normal current flowing through the resistive channel yields a current voltage characteristics 
showing steps at the voltage values given by (11.2). Current driven arrays are somewhat 
more difficult to treat theoretically, because one really has to solve the equation of motion 
(2.4) for a given drive current. Shapiro steps again occur at the voltage values of (11.2). 

Under ^ideal^ conditions, ie in zero temperature and external field, an M x array of 
such junctions should act as a single junction with an applied current I/M. The total voltage 
drop across the array should then be N times the one across each junction and ^giant^ 



Shapiro steps should occur at voltages 

Vn = NVn (11.3) 

Reference 1 gives many an excellent overview over various more elaborate aspects of Shapiro 
steps in J J A, such as the occurrence of fractional giant steps and of suharmonic steps. 

Application for microwave radiation 

A Josephson junction is a voltage controlled oscillator capable of generating electromag- 
netic radiation [122]. The dynamics of the supreconducting phase can be formulated in the 
framework of the capacitively and resistively shunted junction model presented in section 
2.1. A dc bias voltage which excites the damped non-linear oscillator has to be added to the 
equations of motion given there. Quantitative estimates show, however, that the linewidth 
of a single junction is rather large [122]. One obvious way to overcome this drawback is 
now precisely to use one- or two-dimensional arrays of junctions and to have them oscillate 
in phase. In the 2D case treated in this article the rf current flows, for example, along the 
rows of the array, whereby care has to be taken that the phases of the current in each row 
are the same. The line width is then reduced, with respect to a single junction, by the 
number of the array sites. Details about the realisation of such voltage controlled microwave 
oscillators in the frequency range between 100 GHz to 1 THz can be found in reference [122] . 

Many of the concepts developed for arrrays can, of course, also be applied to the de- 
scription of thin superconducting films which are, in some sense, a continuum version of the 
discrete lattice. Thus we have mentioned results found on such films in sections 6 and 8. 
Moreover, high temperature superconductors have more or less strongly anisotropic struc- 
tures of rather weakly coupled lattice planes. Thus, for many purposes the phenomenon of 
superconductivity can be restricted to lattice planes, so that many aspects of two-dimensional 
suprerconductivity described in the present article are again applicable. 
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